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NUMERICAL AND ALGEBRAIC 

Gain in decibels of P2 relative to Pi 

G = 101ogio(^2/Pi). 



To within two percent 
(27r) 

Euler-Mascheroni constant^ 7 = 0.57722 
Gamma Function r(cc + 1) = xT^x): 



^ 2.5; TT^ ^ 10; ?^ 20; 2^° « 10^. 



r(i/6) 


= 5.5663 


r(3/5) 


= 1.4892 


r(i/5) 


= 4.5908 


r(2/3) 


= 1.3541 


r(i/4) 


= 3.6256 


r(3/4) 


= 1.2254 


r(i/3) 


= 2.6789 


r(4/5) 


= 1.1642 


r(2/5) 


= 2.2182 


r(5/6) 


= 1.1288 


r(i/2) 


= 1.7725 = 


r(i) 


= 1.0 



Binomial Theorem (good for | a: |< 1 or o; = positive integer) 

00 

(1 + = ^ Q^x^ = l + ax + 



a(a - 1) 2 , «^(« - - 2) 3 , 
X H X + 



2! 



3! 



fe = 



Rothe-Hagen identity'^ (good for all complex x, y, z except when singular): 

n 

EX /X -\- kz\ y /y + {n — k)z\ 

X + kz ^ k ' y + {n — k)z ^ n — k ' 

k=o 

_ X + y /X + y + nz\ 
X + y + nz ^ n ^ 

Newberger's summation formula^ [good for /x nonintegral, Re (o; + 13) > — 1]: 

00 

( — 1)"' Jq, — ^n{z)J^j^^n{z) 



E 



n= — 00 



n + jjb 



Sm /XTT 



Joc + 'ytJi{z)Jp-^^{z). 
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VECTOR IDENTITIES* 

Notation: /, are scalars; A, B, etc., are vectors; 7" is a tensor; / is the unit 
dyad. 

(1) ABxC = AxBC = BCxA = BxCA = CAxB = CxAB 

(2) A X (B X C) = (C X B) X A = (A ■ C)B - (A ■ B)C 

(3) A X (B X C) + B X (C X A) + C X (A X B) = 

(4) (A X B) ■ (C X D) = (A ■ C)(B ■ D) - (A ■ D)(B ■ C) 

(5) (A X B) X (C X D) = (A X B ■ D)C - (A x B ■ C)D 

(6) V(/^) = V(^/) = fVg + gVf 

(7) V-(/A) = /V-A + A.V/ 

(8) V X (/A) = /V X A + V/ X A 

(9) V-(AxB) = B- VxA-A-VxB 

(10) V x (A x B) = A(V ■ B) - B(V ■ A) + (B ■ V)A - (A ■ V)B 

(11) A X (V X B) = (VB) ■ A - (A • V)B 

(12) V(A • B) = A X (V X B) + B X (V X A) + (A • V)B + (B • V)A 

(13) VV = V- V/ 

(14) V^A = V(V • A) - V X V X A 

(15) V X V/ = 

(16) V ■ V X A = 

If ei, 62, ea are orthonormal unit vectors, a second-order tensor T can be 
written in the dyadic form 

(17) T=^ T,,e,e,- 

In cartesian coordinates the divergence of a tensor is a vector with components 

(18) iV-T)i = J2^idTji/dxj) 

[This definition is required for consistency with Eq. (29)]. In general 

(19) V • (AB) = (V • A)B + (A • V)B 

(20) V-(/T) = V/-T+/V-T 
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Let r = ia; + jy + be the radius vector of magnitude r, from the origin to 
the point x, y, z. Then 

(21) V • r = 3 

(22) V X r = 

(23) Vr = r/r 

(24) V(l/r) = -r/r^ 

(25) V • (r/r^) = 47r5(r) 

(26) Vr = / 

If y is a volume enclosed by a surface S and dS — ndS, where n is the unit 
normal outward from V, 

(27) / dVVf= / dSf 
Jv J s 

(28) / dVV • A = / dS- A 
Jv J s 

(29) / dVV-T= dS T 
Jv J s 

(30) / dVV X A = / dS X A 
Jv J s 

(31) / dV(fV''g-gV^f)= / dS-(fVg-gVf) 
Jv J s 

(32) / dV(A - VxVxB-B-VxVxA) 
J V 

= / dS-(BxVxA-AxVxB) 
J s 

If S is an open surface bounded by the contour C, of which the line element is 
dl, 

(33) / dS X V/ = (h dlf 
J s J c 



(34) / dS-V X A = d) dl- A 

s c 

(35) / (dS X V) X A = (j) dl X A 
J s J c 

(36) / dS ■ (V/ xVg)= d) fdg = - (j) gdf 
J s J c J c 



DIFFERENTIAL OPERATORS IN 
CURVILINEAR COORDINATES^ 



Cylindrical Coordinates 

Divergence 

1 a , , 1 OAs dA^ 

V • A = {rAr) + H 

r dr r dcf) dz 

Gradient 

(V/).-§^; (V/)^ = i|^; (V/), = 

or r ocp oz 



Curl 

(V X A)r = 

(V X A)^ = 



1 dA^ dA 



4> 



r d(f) dz 

dAr dA^ 
dz dr 



1 a , , 1 dAr 
(VxA), = - — (rA<^)---f 
r dr r d<p 



Laplacian 



r dr \ dr J r^ dcf)"^ dz^ 
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Laplacian of a vector 



r o(p r 



Components of (A • V)B 



9r r 9^ dz 



ar r a0 az r 



dB^ As dB^ dB 
(A ■ VB);, = Ar + — ^ + A 



z 



dr ' r d(j) ^ dz 



Divergence of a tensor 

(V • I )r = - — [rTrr) H ^— + 

r or r 0(p oz r 



r or r ocp oz r 



(V . T), = - — (rT.,) + -— ^ + — ^ 
r ar r (70 oz 
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Spherical Coordinates 

Divergence 

V-A= {r^Ar) + ^—- — isindAe) + 



dr r sin 6 86 r sin 6 d(j) 



Gradient 



Curl 



(V/). = 1^; (V/), = i|^; {yf)^ = ^—^ 
or r o& r sm & ocp 



Id 1 dAg 

r sm 6^ at^ r sm & ocp 



(V X A), = — - --^(rA^) 
r sm 6^ ocp r or 



1 d ^ X 1 dAr 

(V xA)^= - —(rAo) -f 

r or r oO 



Laplacian 

^ ~ \ 'd^J ^ r2 sin 9 80 'dO J ^ sin^ 6 

Laplacian of a vector 

(V^A) = A - — - — 5^ - ^ _ 2 



, 2 , 2 2 2 cos 6* SA, 



r-^ 9^ r-^ sin-^ 6 r"^ sin-^ ^ 90 



^ 2 aA^ 2cos0 aA( 



(V^A)^ = V^A^ - -— + -___-L + 



sin^ sin 90 sin^ 9 dcj) 
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Components of (A • V)B 

/A vTrD^ A , Ae dBr , dBr AeEff+A^B, 

(A • VtJjr — -^r + + 



dr r 86 r sin dcj) 



(A ■ VU)g = Ar \ \ \ 

dr r dO r sin d(j) r r 



/A VT^^^ A ^^4> , ^0 A^ dB^ A^Br , cot 6 A^Be 

(A • VB)^ = Ar— 1 — - H "I 1 

or r oU r sin 6 ocp r r 

Divergence of a tensor 

(V- T)r = ^^{r^Trr) + ^—.^{siT^OTer) 
r"^ or r sin oo 



r sin dip 



(V- T)e - — — (r^T,e) + ^-T — (sin^T,,) 
or r sm 6^ at^ 



^ 1 9T<^6i _|_ T'er cot dT^^ 
r sin d(f) r r 



ar r sin ou 



^ 1 dT(p(p ^ T^r _|_ cot 9T(^0 
r sin 9</) r r 
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DIMENSIONS AND UNITS 



To get the value of a quantity in Gaussian units, multiply the value ex- 
pressed in SI units by the conversion factor. Multiples of 3 in the conversion 
factors result from approximating the speed of light c = 2.9979 X 10^° cm/sec 
» 3 X 10^° cm/sec. 



Physical 
Quantity 



Sym- 
bol 



Dimensions 



SI 



Gaussian 



SI 
Units 



Conversion 
Factor 



Gaussian 
Units 



t q 


I 

Tfl L 


72 


Q 


t 


Q 


1/2 


P 


PIH 


tq 


7 

I 




t 


tq' 


1 


mP 


t 


Q 




t 




q 


ml/2 


IH 




m 


m 


P 


P 


Q 


ml/2 


P 


P/'^t 


ml 


ml/2 



Capacitance 
Charge 

Charge 
density 

Conductance 

Conductivity 

Current 

Current 
density 

Density 
Displacement 

Electric field 

Elect ro- 
motance 

Energy 

Energy 
density 



C 



I,i 

JJ 

P 
D 

E 

S, 

Emf 
U, W 

w, e 



Pq 
mP 

~P~ 

m 

1p 



p/^t 

_l/2;l/2 



~P 

m 

Tp 



farad 
coulomb 

coulomb 
/m3 

Siemens 



Siemens 
/m 

ampere 



ampere 
/m2 

kg/m^ 



coulomb 
/m2 

volt / m 

volt 

joule 
joule/m^ 



9 X 10 



11 
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3 X 10'' 
3 X 10^ 
9 X lOii 
9 X 10^ 
3 X 10 
3 X 10^ 
10"^ 

127r X 10^ 



- X 10" 

3 



- X 10" 
3 



10" 
10 



cm 

statcoulomb 

statcoulomb 
/cm3 

cm/sec 
sec~i 

statampere 

statampere 
/cm^ 

g/cm^ 



statcoulomb 
/cm^ 

statvolt / cm 

statvolt 

erg 

erg/cm^ 
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Physical 
Quantity 



Sym- 
bol 



Dimensions 



SI 



Gaussian 



SI 
Units 



Conversion 
Factor 



Gaussian 
Units 



Force 

Frequency 

Impedance 

Inductance 
Length 

Magnetic 
intensity 

Magnetic flux 

Magnetic 
induction 

Magnetic 
moment 

Magnetization 

Magneto- 
motance 

Mass 

Momentum 

Momentum 
density 

Permeability 



F 

Z 

L 
I 

H 

B 

m, n 
M 

M, 
Mmf 

m, M 
P,P 



ml 

1^ 
1 

t 

ml' 
ml' 



It 

ml^ 

~tq 

m 

tq 



1 
It 

1 
t 

m 

ml 
t 

m 

¥t 

ml 



ml 

1^ 
1 

t 
t 

1 



m 



1/2 



t 

1/2 



l^/H 
t 



l^/H 

m 

ml 
t 

m 



newton 

hertz 

ohm 

henry 
meter (m) 

ampere— 
turn/m 

weber 
tesla 

ampere— m^ 

ampere— 
turn/m 

ampere- 
turn 

kilogram 

(kg) 
kg— m / s 

kg/m^— s 
henry /m 



10^ 



- X 10 
9 



11 



11 



- X 10 

9 

10^ 

47r X 10"^ 

10® 
10^ 

10^ 

47r X 10"^ 
47r 

To 

10^ 
10^ 

10-1 



— X 10' 

47r 



dyne 
hertz 
sec/cm 

sec^/cm 

centimeter 
(cm) 

oersted 

maxwell 
gauss 

oersted- 
cm^ 

oersted 

gilbert 

gram (g) 

g— cm / sec 
g/cm^— sec 
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Dimensions 








Physical 
Quantity 


Sym- 
bol 






SI 
Units 


Conversion 
Factor 


Gaussian 
Units 


SI 


Gaussian 


PprTTiittivitv 
Polarization 

Potential 


P 

V,cl> 


mP 
Q 


1 

1/2 

m ' 


fa.ra.fl /m 
coulomb/m^ 

volt 


367r X 10^ 
3 X 10^ 

- X 10 ^ 
3 


statcoulomb 
/cm^ 

statvolt 


mP 


P/^t 
m^l'^pl'^ 


Pn 


f 


Power 


P 


mP 


mP 
P 


watt 


10^ 


erg/sec 


Power 
density 




m 
IP 


m 

Ip 


watt / m^ 


10 


erg/cm^— sec 


J- 1 t;i3i3 HI t; 


rt P 


m 
Jp 


m 


T\Ck C/-' o 1 


1 




Reluctance 
Resistance 


n 

R 


mP 
mP 
tq^ 


1 

y 

t 
1 


ampere— turn 
/weber 

ohm 


47r X 10"^ 

- X 10-^^ 
9 


cm~^ 
sec/cm 


Resistivity 


V,P 


mP 
tq^ 


t 


ohm— m 


- X 10-^ 
9 


sec 


Thermal con- 
Huctivitv 




ml 
~P 


ml 
~P 


watt / m— 
dee- (K) 


10^ 


erg/cm— sec- 
dee (K\ 


Time 

Vector 
potential 


t 

A 


t 

ml 


t 


second (s) 
weber/m 


1 

10^ 


second (sec) 
gauss— cm 


1^ 
I 

4. 
t 


t 

I 

4. 
t 


Velocity 


V 


m/s 


10^ 


cm/sec 


Viscosity 


V, M 


m 
It 


m 
It 


kg/m— s 


10 


poise 


Vorticity 
Work 


C 

w 


1 

t 

mP 

t2 


1 

t 

mP 

t2 


s-^ 
joule 


1 

10^ 


sec~^ 
erg 
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INTERNATIONAL SYSTEM (SI) NOMENCLATURE 



Physical 


Name 


Symbol 


Physical 


Name 


Symbol 


Quantity 


of Unit 


for Unit 


Quantity 


of Unit 


for Unit 


* length 


meter 


m 


electric 


volt 


V 






potential 






*mass 


kilogram 


kg 












electric 


ohm 




*time 


second 


s 


resistance 






current 


ampere 


A 


electric 


Siemens 


s 






conductance 






*temperature 


kelvin 


K 














electric 


farad 


F 


amounu oi 


mole 


mol 


C cl p cl C 1 X cl n C G 






substance 
















magnetic rlux 


weber 


Wb 


*luminous 


candela 


cd 






intensity 






magnetic 


henry 


H 








inductance 






t plane angle 


radian 


rad 












lliclg,llc Lie 


tcoia 


T 

_L 


tsolid angle 


steradian 


sr 


intensity 






frequency 


hertz 


Hz 


luminous flux 


lumen 


Im 


energy 


joule 


J 


illuminance 


lux 


Ix 


force 


newton 


N 


activity (of a 


becquerel 


Bq 








radioactive 






pressure 


pascal 


Pa 


source) 






power 


watt 


W 


absorbed dose 


gray 


Gy 








(of ionizing 






electric charge 


coulomb 


C 


radiation) 







*SI base unit jSupplementary unit 



METRIC PREFIXES 



Multiple 


Prefix 


Symbol 


Multiple 


Prefix 


Symbol 


10-1 


deci 


d 


10 


deca 


da 


10-2 


centi 


c 


102 


hecto 


h 


10-^ 


milli 


m 


10^ 


kilo 


k 


10-^ 


micro 




10^ 


mega 


M 


10-9 


nano 


n 


10^ 


giga 


G 


10-12 


pico 


P 


1012 


tera 


T 


10-15 


femto 


f 


10i5 


peta 


P 


10-18 


atto 


a 


IQis 


exa 


E 
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PHYSICAL CONSTANTS (SI)^ 



Physical Quantity 


Symbol 


Value 


Units 


xjOiLzmann consLanx 


I- 

K 


1.3807 


X 


10" 


23 


J IV 


Elementary charge 


e 


1.6022 


X 


10" 


19 




H/iection mass 


me 


9.1094 


X 


10" 


31 


kg 


Proton mass 


mp 


1.6726 


X 


10" 


27 


kg 


Gravitational constant 


G 


6.6726 


X 


10" 


11 


m"^s kg ^ 


Planck constant 


h 


6.6261 


X 


10" 


34 


Js 




h = h/2'K 


1.0546 


X 


10" 


34 


Js 


Speed of light in vacuum 


c 


2.9979 


X 


10* 




m s""*^ 


Permittivity of 




8.8542 


X 


10" 


12 


Fm"^ 


free space 














Permeability of 


MO 


47r 


X 


10" 


7 


H m" ^ 


free space 














Proton/electron mass 


mp/me 


1.8362 


X 


10^ 






ratio 














Electron charge/mass 


e/nie 


1.7588 


X 


10^1 


Ckg"i 


ratio 














Rydberg constant 


me 
oeo cn,-^ 


1 new A 


X 


10^ 




m"^ 


Bohr radius 


u2 1 2 

ao = eoAi /Tzme 


5.2918 


X 


10" 


11 


m 


Atomic cross section 


2 

Trao 


8.7974 


X 


10" 


21 


m^ 


Classical electron radius 


2 / . 2 

Te = e /47reomc 


2.8179 


X 


10" 


15 


m 


Thomson cross section 


(87r/3)re^ 


6.6525 


X 


10" 


29 


m^ 


Compton wavelength of 


h/rrieC 


2.4263 


X 


10" 


12 


m 


electron 


h/rrieC 


3.8616 


X 


10" 


13 


m 


Fine-structure constant 


q; = e /2eohc 


7.2972 


X 


10" 


3 








137.038 






First radiation constant 


Ci = 27rhc 


3.7418 X 


10" 


16 


Wm^ 


Second radiation 


C2 = hc/fc 


1.4388 X 


10" 


2 


mK 


constant 














Stefan-Boltzmann 


cr 


5.6705 X 


10" 


8 


Wm"2K"^ 


constant 
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Physical Quantity 


Symbol 


Value 


Units 


Wavelength associated 


Ao = 


- hc/e 


1.2398 X 10" 


6 


m 


W 1 Lli -L t; V 












Frequency associated 




: e/h 


2.4180 X 10^^ 


XT 

Hz 


W i t li -L V 












Wave number associated 


ko - 


-ejhc 


8.0655 X 10^ 




m 


1 1", Vl 1 P \l 

VV J.UJ.± _l_ V y V 












Energy associated with 


hh>o 




1.6022 X 10" 


19 


J 


1 pV 

J. c V 












Energy associated with 


he 




1.9864 X 10" 


25 


T 
J 


1 m ~ ■'^ 
J. ill 












-L;11C1 dooLiC/lclLt^d W 1 til 


Q 

me 


/8eo 


13.606 




pV 

c V 


1 Rydberg 












Energy associated with 


k/e 




8.6174 X 10" 


5 


eV 


1 Kelvin 












Temperature associated 


e/ k 




1.1604 X 10^ 




K 


with 1 eV 












Avogadro number 


Na 




6.0221 X 10^^ 


mol""*^ 


Faraday constant 




Map 


9.6485 X 10^ 




Cmol"^ 


Gas constant 


K = 


NAk 


8.3145 




JK" Vol"^ 


Loschmidt's number 


no 




2.6868 X 10^^ 


m"^ 


(no. density at STP) 












Atomic mass unit 






1.6605 X 10" 


27 


kg 


Standard temperature 


-t 




273.15 




K 


Atmospheric pressure 


Po — 


- llQhul 


1.0133 X 10^ 




Pa 


Pressure of 1 mm Hg 






1.3332 X 10^ 




Pa 


(1 torr) 












Molar volume at STP 


Vo = 


= RTq/po 


2.2414 X 10" 


2 


m^ 


Molar weight of air 






2.8971 X 10" 


2 


kg 


calorie (cal) 






4.1868 




J 


Gravitational 


9 




9.8067 




m s"'^ 


acceleration 













15 



PHYSICAL CONSTANTS (cgs)^ 



Physical Quantity 


Symbol 


Value 


Units 


Boltzmann constant 


k 


1.3807 


X 


10" 


16 


erg/ deg (K) 


Elementary charge 


e 


4.8032 


X 


10" 


10 


statcoulomb 














(statcoul) 


Electron mass 


me 


9.1094 


X 


10" 


28 


g 


Proton mass 


nip 


1.6726 


X 


10" 


24 


g 


Gravitational constant 


G 


6.6726 


X 


10" 


8 


dyne-cm^ /g'^ 


Planck constant 


h 


6.6261 


X 


10" 


27 


erg-sec 




h = h/2'K 


1 0546 


V 


10" 


27 


erg-sec 


Speed of light in vacuum 


c 


2 9979 


V 
/\ 


10^0 


cm/sec 


Proton/electron mass 


mp/rue 


1 8362 




10^ 






ratio 














Electron charge/mass 


e/rue 


5.2728 


X 


10^^ 


statcoul/g 


ratio 














Rydberg constant 


27r^me^ 












~ ch^ 


1.0974 


X 


10^ 




cm" ^ 


Bohr radius 


ao = n /me 


5.2918 


X 


10" 


9 


cm 


Atomic cross section 




8.7974 


X 


10" 


17 


cm^ 


Classical electron radius 


1 1 2 

Te = e /mc 


2.8179 


X 


10" 


13 


cm 


Thomson cross section 


(87r/3)re^ 


6.6525 


X 


10" 


25 


cm^ 


Compton wavelength of 


h/m^c 


2.4263 


X 


10" 


10 


cm 


electron 


h/me,c 


3.8616 


X 


10" 


11 


cm 


Fine-structure constant 


a = /he 


7.2972 


X 


10" 


3 






a-' 


137.038 






First radiation constant 


Ci = 27rhc^ 


3.7418 X 


10" 


5 


erg-cm^ /sec 


Second radiation 


C2 = hc/k 


1.4388 




cm-deg (K) 


constant 














Stefan-Boltzmann 


a 


5.6705 X 


10" 


5 


erg/cm^- 


constant 












sec-deg'* 


Wavelength associated 


Ao 


1.2398 X 


10" 


4 


cm 


with 1 eV 
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Physical Quantity 


Symbol 


Value 


Units 


Frequency associated 






2.4180 X 


10^4 


Hz 


W 1 tli -L C V 












Wave number associated 


ko 




8.0655 X 


10^ 


cm 


wif.Vi 1 pV 

VV lljXJ. -L V 












Energy associated with 






1.6022 X 


10-12 


erg 


1 pV 

J. c V 












Energy associated with 






1.9864 X 




erg 


1 cm 


















13.606 


pV 

t; V 


1 Rydberg 












Energy associated with 






8.6174 X 


10-^ 


eV 


1 deg Kelvin 












Temperature associated 






1.1604 X 


10^ 


deg (K) 


with 1 eV 












Avogadro number 


Na 




6.0221 X 


10^3 


mol-i 


Faraday constant 


/ — 


Map 


2.8925 X 


10^4 


statcoul/mol 


Gas constant 


H = 


AT 


8.3145 X 


10^ 


erg/ deg-mol 


Loschmidt's number 


no 




2.6868 X 


10^9 


cm — ^ 


(no. density at STP) 












Atomic mass unit 


Tn„, 

1 1 On, 




1.6605 X 


10-24 


g 


Standard temperature 


-L 




273.15 


deg (K) 


Atmospheric pressure 


Po — 


ItQtvl 


1.0133 X 


10*^ 


dyne/cm2 


Pressure of 1 mm Hg 






1.3332 X 


10^ 


dyne/cm2 


(1 torr) 












Molar volume at STP 


Vo = 


RTo/po 


2.2414 X 


10^ 


cm^ 


Molar weight of air 


Mair 




28.971 


g 


calorie (cal) 






4.1868 X 


10" 


erg 


Gravitational 


9 




980.67 


cm/sec2 


acceleration 
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FORMULA CONVERSION^ 

Here a = 10^cmm"\ f3 = 10^ergJ"\ eo = 8.8542 X 10"^^Fm"\ 
fXQ = 47rXlO"^Hm-\ c= (eoA^o)"^/^ = 2.9979 X 10® m s" ^ and h = 1.0546X 
10 J s. To derive a dimensionally correct SI formula from one expressed in 
Gaussian units, substitute for each quantity according to Q = kQ, where k is 
the coefficient in the second column of the table corresponding to Q (overbars 

— 2 2 

denote variables expressed in Gaussian units). Thus, the formula ao = h /me 
for the Bohr radius becomes aao = (?j,/3)^/[(m/3/a^)(e^CK/3/47reo)], or ao = 
eoh^ /TTme^ . To go from SI to natural units in which ^ = c = 1 (distinguished 
by a circumflex), use Q — k~^Q, where k is the coefficient corresponding to 
Q in the third column. Thus do = 4:7reoh^ /[(rhh/c)(e'^eohc)] = 47r/me^. (In 
transforming from SI units, do not substitute for eo, fj^o, or c.) 



Physical Quantity 


Gaussian Units to SI 


Natural Units to SI 


Capacitance 


CK / 47reo 




Charge 


(a/3/47ren)^/^ 


(eoM"'/^ 


Charge density 


(/3/47ra^eo)^/^ 


(eoM"'/^ 


Current 


(a/3/47reo)^/^ 


\ t W / / 


Current densitv 






Electric field 


(47r^eo/a3)i/2 


(eo/M'/' 


Electric potential 


(47r/3eo/a)i/2 


(eoAc)i/2 
60-' 


Electric conductivity 


(47reo)-' 


Energy 


/3 




Energy density 






Force 


/3/a 


(he)-' 


Frequency 


1 


Inductance 


47reo / Oi 




Length 


a 


1 


Magnetic induction 


(47r/3/aVo)'/^ 




Magnetic intensity 


(47r/xo/3/a=^)^/2 


(/xoAc)i/2 


Mass 


c/^i 


Momentum 


P/a 




Power 






Pressure 




(he)-' 


Resistance 


47reo / Oi 


(6o/Mo)'/" 


Time 


1 


c 


Velocity 


a 
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MAXWELL'S EQUATIONS 



Name or Description 



SI 



Gaussian 



Faraday's law 



Ampere's law 

Poisson equation 

[Absence of magnetic 
monopoles] 

Lorentz force on 
charge q 

Constitutive 
relations 



V X E = -■ 



as 



V X H = + J 

dt 

V • D = p 

V • B = 

g (E + V X B) 

D = eE 
B = /iH 



1 aB 

V X E = 

c dt 

1 aD 47r 

V X H = + — J 

c dt c 

V • D = 47rp 

V • B = 



g ( E + -V X B 



D = eE 
B = /xH 



In a plasma, jjL jjLo = 47r X 10 Hm (Gaussian units: jjL ~ 1). The 

permittivity satisfies e ?a eo = 8.8542 x 10~^^Fm~^ (Gaussian: e ?a 1) 
provided that all charge is regarded as free. Using the drift approximation 
Vj^ = E x B/S'^ to calculate polarization charge density gives rise to a dielec- 
tric constant K = e/eo = l + 367r X l(fp/B^ (SI) = l + 47rpc^/S^ (Gaussian), 
where p is the mass density. 

The electromagnetic energy in volume V is given by 



W = - I dV{H • B + E • D) 



(SI) 



= — / dV(n B + E D) 
87r 



(Gaussian). 



Poynting's theorem is 



dW 
~dt 



+ 



N ■ dS = - 



dVJ ■ E, 



where S is the closed surface bounding V and the Poynting vector (energy fiux 
across S) is given by N = E x H (SI) or N = cE x H/47r (Gaussian). 
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ELECTRICITY AND MAGNETISM 

In the following, e = dielectric permittivity, /u, = permeability of conduc- 
tor, fx' = permeability of surrounding medium, cr = conductivity, / = ui /27r = 
radiation frequency. Km = m/mo and Ke = e/eo- Where subscripts are used, 
'1' denotes a conducting medium and '2' a propagating (lossless dielectric) 
medium. All units are SI unless otherwise specified. 

Permittivity of free space eo = 8.8542 x 10~^^ Fm~^ 

Permeability of free space fxo = 47r X 10-"^ H m 

= 1.2566 X 10"^ Hm"^ 

Resistance of free space Rq = (/xo/eo)^/^ = 376.73 Q 

Capacity of parallel plates of area C = eA/d 

A, separated by distance d 

Capacity of concentric cylinders C = 27reZ/ ln(6/a) 

of length I, radii a, b 

Capacity of concentric spheres of C — 4'Keab/ (b — a) 
radii a, b 

Self-inductance of wire of length L = ij,1/87t 
I, carrying uniform current 

Mutual inductance of parallel wires L = (/x'Z/47r) [1 + 41n(d/a)] 
of length I, radius a, separated 

by distance d 

Inductance of circular loop of radius L — b ^/j,' [ln(86/a) — 2] + A''/4j' 
b, made of wire of radius a, 
carrying uniform current 

Relaxation time in a lossy medium r = e/cr 

Skin depth in a lossy medium S = (2/c<j/icr)"'"/^ = (tt/z/ct)""'"/^ 

Wave impedance in a lossy medium Z = + icr /to)]'^^'^ 

Transmission coefficient at T = 4.22 x lQ~^{f Kmii^e^/ cy)^^^ 

conducting surface^ 
(good only for T <^ 1) 

Field at distance r from straight wire Be — yLiI/27rr tesla 

carrying current / (amperes) = 0.2J/r gauss (r in cm) 

Field at distance z along axis from Bz = //a^//[2(a^ + z^)"^/^] 

circular loop of radius a 
carrying current / 
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ELECTROMAGNETIC FREQUENCY/ 
WAVELENGTH BANDS^° 



Designation 


Frequency Range 


Wa,vplpnef 


t li R a,n ffp 


Lower 


Upper 


Lower 


Upper 


ULF* 




30 Hz 


10 Mm 




VF* 


30 Hz 


300 Hz 


1 Mm 


10 Mm 


ELF 


300 Hz 


3 kHz 


100 km 


IMm 


VLF 


3 kHz 


30 kHz 


10 km 


100 km 


LF 


30 kHz 


300 kHz 


1km 


10 km 


MF 


300 kHz 


3 MHz 


100 m 


1km 


HF 


3 MHz 


30 MHz 


10 m 


100 m 


VHF 


30 MHz 


300 MHz 


Im 


10 m 


UHF 


300 MHz 


3 GHz 


10 cm 


Im 


SHFf 


3 GHz 


30 GHz 


1 cm 


10 cm 


S 


2.6 


3.95 


7.6 


11.5 


G 


3.95 


5.85 


5.1 


7.6 


J 


5.3 


8.2 


3.7 


5.7 


H 


7.05 


10.0 


3.0 


4.25 


X 


8.2 


12.4 


2.4 


3.7 


M 


10.0 


15.0 


2.0 


3.0 


P 


12.4 


18.0 


1.67 


2.4 


K 


18.0 


26.5 


1.1 


1.67 


R 


ZD.O 


/inn 


0.75 


1.1 


EHF 


30 GHz 


300 GHz 


1 mm 


1 cm 


Submillimeter 


300 GHz 


3 THz 


100 /im 


1 mm 


Infrared 


3THz 


430 THz 


700 nm 


100 /xm 


Visible 


430 THz 


750 THz 


400 nm 


700 nm 


Ultraviolet 


750 THz 


30PHz 


10 nm 


400 nm 


X Ray 


30PHz 


3EHz 


100 pm 


10 nm 


Gamma Ray 


3EHz 






100 pm 



In spectroscopy the angstrom is sometimes used (lA = 10 ^ cm = 0.1 nm). 
*The boundary between ULF and VF (voice frequencies) is variously defined. 
fThe SHF (microwave) band is further subdivided approximately as shown.^^ 
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AC CIRCUITS 

For a resistance R, inductance L, and capacitance C in series with 
a voltage source V = Vq exp(iu;t) (here i = the current is given 

by J = dq/dt, where q satisfies 

d^q dq q 
L—^ + R— + — = V. 
dt^ dt C 

Solutions are q{t) = Qa + Qt, I{t) = Is + It, where the steady state is 
Is = iujqs = V/Z in terms of the impedance Z = R + i(ujL — 1/ujC) and 
It = dqt/dt. For initial conditions ^(O) = go = + Qs, 1(0) = -^o, the 
transients can be of three types, depending on A = R^ — 4L/C: 

(a) Overdamped, A > 



Qt 
It 



Jo + 7+ go 



7+ - 7- 
7+ (Jo + 7-go) 
7+ - 7- 



exp(-7_t) 



exp(-7+t) 



Jo +7-^0 . , . 
exp(-7+t), 

7+ - 7- 

7_(Jo + 7+go) 



7+ - 7- 



exp(-7_t), 



where j± = (R ± A^/^)/2L; 
(b) Critically damped, A = 



Qt = [qo + (lo + lRQo)t] exp(-7Ht), 
It = [lo - (lo + 7hQo)7h*] exp(-7Ht), 



where 7h — R/2L; 

(c) Underdamped, A < 



If = 



IRQO + Iq 



sin (jjxt -\- qo cos uxt 



exp(-7Ht), 



iQ cosujit sm(a;ir) 



UJl 



exp(-7Rt), 



Here cui = uJo{l — J?^C/4L)''"/^, where cuq = (LC)""*"/^ is the resonant 
frequency. At uj — ujq, Z = R. The quality of the circuit is Q — ujqL/R. 
Instability results when L, R, C are not all of the same sign. 
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DIMENSIONLESS NUMBERS OF FLUID MECHANICS^^ 



Naine(s) 


Symbol 


Definition 


S ignificance 


Alfven, 


Al, Ka 


Va/V 


* (Magnetic force/ 


Karman 






, . T r \ 1 /2 

mertial force) ' 


Bond 


Bd 


{p' - p)L^g/T. 


Gravitational force/ 






surface tension 


Boussinesq 


B 




(Inertial force/ 

\ 1 / 2 

gravitational force) ^ 


Brinkman 


Br 




Viscous heat/conducted heat 


Capillary 


Cp 




Viscous force/surface tension 


Carnot 


Ca 


{T2-T^)/T'2 


Theoretical Carnot cycle 








Cauchy, 


Cy, Hk 


pV^ /T - 


Inertial force/ 


rlooke 




compressibility force 


Chandra- 


Ch 


/ pi/f] 


Magnetic force/dissipative 


sekhar 




forces 


v^iaUblUo 




LV^p/kAT 


xvinetic eneigy now rare/ neat 
conduction rate 


Cowling 


C 


(VA/Vf =Al^ 


Magnetic force/inertial force 


Crispation 


v^r 


,, u~ /V T 
fJLK / Z-iL/ 


Effect of diffusion/effect of 
surface tension 


Dean 


D 




Transverse flow due to 






curvature/ longitudinal flow 


[Drag 


Cd 


{p' — p)Lg/ 
pW 


Drag force/inertial force 


coefficient] 






Eckert 


E 


V'^/cpAT 


Kinetic energy/change in 






thermal energy 


Ekman 


Ek 


(Ro/Re)^/^ 


(Viscous force/ Coriolis force) 


Euler 


Eu 


Ap/pV^ 


Pressure drop due to friction/ 






dynamic pressure 


Froude 


Fr 


-1 / o 

V/igL)'/^ 
V/NL 


1 (Inertial force/gravitational or 
buoyancy force) ''"/'^ 


Gay— Lussac 


Ga 


1//3AT 


Inverse of relative change in 




volume during heating 


Grashof 


Gr 


gL^/3AT/u^ 


Buoyancy force/ viscous force 


[Hall 


Ch 




Gyrofrequency / 


coefficient] 






collision frequency 



*(t) Also defined as the inverse (square) of the quantity shown. 
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Name(s) 


Symbol 


Definition 


Significance 


Hartmann 


H 


BL/it^T])'/^ = 
(Rm ReC)^/^ 


(Magnetic force/ 
dissipative force) 


Knudsen 


Kn 


X/L 


Hydrodynamic time/ 
collision time 


Lewis 


Le 


k/V 


*Thermal conduction/molecular 

diffusion 


Lorentz 


Lo 


V/c 


Magnitude of relativistic effects 


Tiiinrlnni sf, 


Lu 


U n Tj \/ A Ti — 

Al Rm 


-T y R forpp /rpsistivp m^iPTiptip 
diffusion force 


Mach 


M 


V/Cs 


Magnitude of compressibility 

effects 


Magnetic 
Mach 


Mm 


V/Va = A\-^ 


(Inertial force/magnetic force) "'^^'^ 


Magnetic 
R pvnolHs 


Rm 


IJ.oLV/ri 


Flow velocity/magnetic diffusion 

VP 1 o r i t V 

V V J. V y V - ± U y 


Newton 


Nt 


F/pL'^V'^ 


Imposed force/inertial force 


Nusselt 


N 


otL/k 


Total heat transfer/thermal 

C011C1U.CL1011 


Peclet 


Pe 


LV/i<i 


Heat convection/heat conduction 


Poisseuille 


Po 




Pressure force/viscous force 


Prandtl 


Pr 


V 1 K 


Momentum diffusion/ 
iieaij Qinubioii 


Rayleigh 


Ra 




Buoyancy force/difTusion force 


Reynolds 


Re 


LV/v 


Inertial force/ viscous force 


Richardson 


Ri 


(NH/AV)^ 


Buoyancy effects/ 

vertical shear effects 


Rossby 


Ro 


y/2QLsin A 


Inertial force/Cor iolis force 


Schmidt 


Sc 


u/V 


Momentum diffusion/ 
molecular difTusion 


Stanton 


St 


/ T 7" 

a/pCpV 


ihermai conduction loss/ 
heat capacity 


Stefan 


Sf 


crLT^/k 


Radiated heat/conducted heat 


Stokes 


s 




Viscous damping rate/ 
vibration frequency 


Strouhal 


Sr 


fL/V 


Vibration speed/flow velocity 


Taylor 


Ta 


.<Sl/v) 


Centrifugal force/viscous force 
r Oentrifueral force/ 
viscous force) ""^/^ 


Thring, 
Boltzmann 


Th, Bo 


pCpV/eaT^ 


Convective heat transport/ 
radiative heat transport 


Weber 


W 


pLV'^/H 


Inertial force/surface tension 
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Nomenclature: 

B Magnetic induction 

Cs,c Speeds of sound, light 

Cp Specific heat at constant pressure (units m.^ K~^) 

D = 2R Pipe diameter 

F Imposed force 

/ Vibration frequency 

g Gravitational acceleration 

H, L Vertical, horizontal length scales 

k = pCpK Thermal conductivity (units kgm~^ 

N = (g / H)'^^'^ Brunt— Vaisala frequency 

R Radius of pipe or channel 

r Radius of curvature of pipe or channel 

vl Larmor radius 

T Temperature 

V Characteristic flow velocity 

Va = S/(/xop)^/^ Alfven speed 

dT 

a Newton's-law heat coefficient, k—— = a AT 

ox 

(3 Volumetric expansion coefficient, dV/V — /3dT 

r Bulk modulus (units kgm~^ 

AR, AV, Ap, AT Imposed differences in two radii, velocities, 

pressures, or temperatures 

e Surface emissivity 

rj Electrical resistivity 

K, I> Thermal, molecular diffusivities (units m^s~^) 

A Latitude of point on earth's surface 

A Collisional mean free path 

jjL = pv Viscosity 

jjio Permeability of free space 

V Kinematic viscosity (units m^s~^) 
p Mass density of fluid medium 

p' Mass density of bubble, droplet, or moving object 

E Surface tension (units kgs~^) 

cr Stefan— Boltzmann constant 

Q. Solid-body rotational angular velocity 
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SHOCKS 



At a shock front propagating in a magnetized fluid at an angle with 
respect to the magnetic induction B, the jump conditions are ^^'^^ 

(1) pU = pU = q; 

(2) pU''+p + B^/2p, = pU^+p + B^/2p,; 

(3) pUV - B\\B^/p = pUV - B\\B^/p- 

(4) ^11 

(5) UB^ - VB\\ = UB^ - 

(6) ^(U^ + V'')+w+iUBj^ - VB\\B^)/l^pU 

= + V^) + w + {UBj^ - VB\\B^)/iJipU. 

Here U and V are components of the fluid velocity normal and tangential to 
the front in the shock frame; p = 1/v is the mass density; p is the pressure; 
B±_ = Ssin^, B\\ = Bcos0\ p is the magnetic permeability {p = 47r in cgs 
units); and the speciflc enthalpy is tt; = e + pv, where the speciflc internal 
energy e satisfles de = Tds — pdv in terms of the temperature T and the 
speciflc entropy s. Quantities in the region behind (downstream from) the 
front are distinguished by a bar. If B = 0, then^^ 

(7) U -U =[{p-p){v-v)]'/^; 

(8) (p-p)(^-^)-i = g2. 

(9) w — w = \{p — p){v + v); 

(10) e — e = ^(p + p){v — v). 

In what follows we assume that the fluid is a perfect gas with adiabatic index 
7 = 1 + 2/n, where n is the number of degrees of freedom. Then p = pRT/m^ 
where R is the universal gas constant and m is the molar weight; the sound 

speed is given by Cs"^ = {dp/dp)s = ^pv; and w = 76 = 7Pf/(7 — 1)- For a 

general oblique shock in a perfect gas the quantity X — (U/Va)^ satisfles''^'* 

(11) (X - 13/ a) (X-cos'^ = Xsin^ {[1 + (r - l)/2a]X - cos^ 6^ where 
r = p/p, a = i [7 + 1 - (7 - l)r], and /3 = C^V^A^ = 4:7vyp/B^. 

The density ratio is bounded by 

(12) Kr < (^+l)/(^-l). 

If the shock is normal to B (i.e., if ^ = 7r/2), then 

(13) = (r/a) {Cs^ + Va^ [1 + (1 - 7/2)(r - 1)]}; 

(14) U/U = B/B = r; 
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(15) V = V; 

(16) p = p + (1 - r~^)pU'^ + (1 - r^)B'^/2fx. 

If ^ = 0, there are two possibihties: switch-on shocks, which require /3 < 1 and 
for which 

(17) = tVa'^; 

(18) U = Va^/U; 

(19) S^2 = 2S||2(r-l)(a-/3); 

(20) y = f7Sx/S||; 

(21) p = p + pU^{l - a + /3)(1 - r"^), 

and acoustic (hydrodynamic) shocks, for which 

(22) C/2 = (r/a)C3^ 

(23) U = L//r; 

(24) y = 5x = 0; 

(25) p = p + pU'^il - r-i). 

For acoustic shocks the specific volume and pressure are related by 

(26) v/v = [(7 + l)p + (7 - l)p] I [(7 - \)p + (7 + \)p\. 
In terms of the upstream Mach number M = U /C a-, 

(27) pIp = v/v = U/U = (7 + 1)mV[(7 - 1)^^ + 2]; 

(28) p/p = (27M2 - 7 + l)/(7 + 1); 

(29) T/T = [(7 - 1)M2 + 2](27M2 - 7 + l)/(7 + ifM'^; 

(30) = [(7 - 1)M2 + 2]/[27M2 -7 + 1]- 
The entropy change across the shock is 

(31) As = s-s = c^ ln[ip/p)ip/py], 

where c^, = R/ (7 — l)m is the specific heat at constant volume; here R is the 
gas constant. In the weak-shock limit (M — >■ 1), 

/QON A ^ 27(7 - 1) ^,^2 in3_ 3 

(32) As ^ c^————(M -1) — _-_(M - 1) . 

3(7 + 1) 3(7 + l)m 

The radius at time t of a strong spherical blast wave resulting from the explo- 
sive release of energy E in a medium with uniform density p is 

(33) Rs = Co{Et^/pV/^ 

where Co is a constant depending on 7. For 7 — 7/5, Co = 1.033. 
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FUNDAMENTAL PLASMA PARAMETERS 



All quantities are in Gaussian cgs units except temperature (T, Tg, Tj) 
expressed in eV and ion mass (m^) expressed in units of the proton mass, 
fj, = mi/mp-, Z is charge state; k is Boltzmann's constant; K is wavenumber; 
7 is the adiabatic index; In A is the Coulomb logarithm. 

Frequencies 

electron gyrofrequency 



ion gyrofrequency 
electron plasma frequency 

ion plasma frequency 



electron trapping rate 

ion trapping rate 

electron collision rate 
ion collision rate 
Lengths 

electron deBroglie length 

classical distance of 
minimum approach 

electron gyroradius 
ion gyroradius 

electron inertial length 
ion inertial length 
Debye length 



fpe — 



UJpe — 



fpi — 



fee = UJcef^TT = 2.80 X 10^ B Hz 

ujce = eB/rueC = 1.76 X lO^Srad/sec 
fci = UJai/'iTT = 1.52 X lO^Zfx'^Bllz 

(jjci = ZeB/ruiC = 9.58 X 10^Z/i~^S rad/sec 
U}pe/27T = 8.98 X lO^ne^/^ Hz 
(47rneeVrne)'/^ 
5.64 X lO^ne^/^ rad/sec 
u;pi/27r 

2.10 X lO^Z/x-^/^ni^/^Hz 

1.32 X 10^Z;u"^/^ni^/^rad/sec 

7.26 X lO^K^/^^^/^ sec"^ 
(ZeKE/rmy/^ 
= 1.69 X 10^zV2Ki/2^i/2^-i/2g^^-i 

i/p =2.91 X 10"^np In ATp"^/^ 



00 



sec 



Ui =4.80 X lO-^Z^/x-^/^nilnATi-^/^sec"^ 

A = h/{mekTey^'^ = 2.76 X 10"^Te"^/^cm 
e^/fcT = 1.44 X 10"^T-^ cm 

re = VTe/i^ce = 2 .38Te^ B' ^ cm 

Ti — VTi/uJci 

= 1.02 X 10^/i^/^Z"^T^^/^S"^ cm 



c/u> 



pe 



5.31 X lO^n -^/^ c 



m 



c/ujpi = 2.28 X 10^Z"^(/x/ni)^/^ cm 

= (fcT/47rne^)^/2 = 7.43 x lO^T^/^n"^/^ 
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Velocities 

electron thermal velocity 

ion thermal velocity 
ion sound velocity 

Alfven velocity 

Dimensionless 

(electron/proton mass ratio) ""^^^ 

number of particles in 
Debye sphere 

Alfven velocity /speed of light 

electron plasma/gyrofrequency 
ratio 

ion plasma/gyrofrequency ratio 
thermal/magnetic energy ratio 
magnetic/ion rest energy ratio 

Miscellaneous 

Bohm diffusion coefficient 

transverse Spitzer resistivity 



VTe = (kTe/me)^^^ 

= 4.19 X lO^Te^/^ cm/sec 
VTi = {kTi/rmy/^ 

= 9.79 X lO^/i-^/^Ti^/^ cm/sec 
Cs = {nZhT^lruifl'^ 

= 9.79 X 10^ (7ZTe///)^/^ cm/sec 
VA = B/(4'Knimiy/'^ 

= 2.18 X 10^ V"^^^^i"^'^^-B cm/sec 

(me/mp)^/^ = 2.33 x 10"^ = 1/42.9 
(47r/3)nAo^ = 1.72 x lO^T^/^n"^/^ 



va/c = 7.28Ai"^/^ni"^/^S 



3.21 X 10"^ne^/^S" 



/3 = SnnkT/B'^ = 4.03 x IQ-^^nTB-^ 



1^2 



Db = (ckT/16eB) 

= 6.25 X lO^TS"^ cm^/sec 



r]± = 1.15 X 10"^^ZlnAT" "'"sec 
= 1.03 X 10"^Zln AT"^/^ Qcm 



-3/2 



The anomalous collision rate due to low-frequency ion-sound turbulence is 

u* » ojpeW/kT - 5.64 X lO^ne^^^M^/fcT sec"\ 

where W is the total energy of waves with uj/K < VTi- 
Magnetic pressure is given by 



mag 



= B'^/Stt = 3.98 X 10^(B/Bof dynes/cm^ = 3.93{B/Bof atm, 



where = lOkG = 1 T. 

Detonation energy of 1 kiloton of high explosive is 



Wi^T = 10^^ cal = 4.2 X 10^^ erg. 
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PLASMA DISPERSION FUNCTION 

Definition^^ (first form valid only for Im ^ > 0): 

, r+°° dt exp f-t^) / X /"^ /X 



- oo 



Physically, C, = x + iy is the ratio of wave phase velocity to thermal velocity. 
Differential equation: 

dZ , X . X 1/2 d^Z dZ 

= -2 (1 + CZ) , Z(0) - iTT^/^; + 2C + 2Z = 0. 

Real argument (y = 0): 

Z(cc) = exp {^—x^^ ^^Tz^^^ — 2 dt exp 

Imaginary argument {x = 0): 

Z{iy) = iTT^/^ exp (y^^ [1 - erf(y)] . 
Power series (small argument): 

Z(C) = iTT^/^ exp (-C^) - 2C (l - 2C^/3 + 4C^/15 - 8C®/105 + ■ ■ ■) 
Asymptotic series, |C| ^ 1 (Ref. 17): 

Z(C) = iTT^/'a exp (-C") - C"' (l + 1/2C' + 3/4C^ + 15/8C^ + • • •) 
where 





y 




{■ 


\y\ 


< \x\ 




y 


< -[■ 



Symmetry properties (the asterisk denotes complex conjugation): 

Z(C*) = -[Z(-C)]*; 

Z(C*) = [Z{C)] * + 2z7r^/" exp[-(C*)'] (2/ > 0). 

Two-pole approximations^® (good for ^ in upper half plane except when y < 
TT^/^cc^ exp(-a;^), x > 1): 

, , 0.50 + 0.81i 0.50 - O.Sli 
^ C ~ 7 a = 0.51-0.8H; 

a - C a* + C 

, , , 0.50 + 0.96i 0.50 - 0.96i 

Z\0^—, h rr^, 6 = 0.48 - 0.91i. 

(6-C)2 (?>* + C)^ 
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COLLISIONS AND TRANSPORT 



Temperatures are in eV; the corresponding value of Boltzmann's constant 
is k = 1.60 X lO"''"'^ erg/eV; masses yu, /u' are in units of the proton mass; 
Coi = Za^e is the charge of species a. All other units are cgs except where 
noted. 

Relaxation Rates 

Rates are associated with four relaxation processes arising from the in- 
teraction of test particles (labeled a) streaming with velocity Va through a 
background of field particles (labeled /3): 



slowing down 



transverse diffusion 



parallel diffusion 



energy loss 



dt 

(Vc 

dV 

d . 



s 



- x2 a\/3 2 



- x2 oc\l3 2 

— Vq,)|| = i^ii Va 



at 



where Vq, = |vq,| and the averages are performed over an ensemble of test 
particles and a Maxwellian field particle distribution. The exact formulas may 
be written"*^^ 



"'^'^ — (1 + moc/mi3)ijj(x°'^'^)u^^'^ ; 

2 [(1 - l/2a;"\^)V^(a;"\^) + V^'(a;"\^)] u^^^ ■ 



«\/3 



oc\l3 



oc\l3 



= 2 



a\/3 
' 



where 



a\B , / Q:\P -^/r>;rn 



2 _ 2 X 2.. 3. _,Q:\/3 2 



V^(a;) = / dtt^^^e *; V^'(a;) = -^, 

and Xctj3 = IuAq,^ is the Coulomb logarithm (see below). Limiting forms of 
Us, u± and un are given in the following table. All the expressions shown 



31 



have units cm^sec~^. Test particle energy e and field particle temperature T 
are both in eV; = mi/rrip where rrip is the proton mass; Z is ion charge 
state; in electron— electron and ion— ion encounters, field particle quantities are 
distinguished by a prime. The two expressions given below for each rate hold 

for very slow {x°^'^^ <C 1) and very fast (a:"^^ ^ 1) test particles, respectively. 

Slow Fast 

Elect ron— elect ron 



e 1 e / \ 




5.8 X 


10 


-6^-3/2 


7.7 


X 


10" 


6 -3/2 

e 






5.8 X 


10 


-6^-1/2^-1 


7.7 


X 


10" 


6 -3/2 


I^" 1 / TleXee 




2.9 X 


10 


-6^,-1/2^-1 


3.9 


X 


10" 


6rp^-5/2 


Electron— ion 


















e\i / ryl \ 

I/g' lUiZ Xei 




0.23/x 


3/2 


j,-3/2 


3.9 


X 


10" 


6 -3/2 


e\i / ryl \ 

V ^ /riiZ Xei 


rsj 


2.5 X 


10 


-4^1/2^-1/2^-1 


7.7 


X 


10" 


6^-3/2 


i^ll IriiZ Xei 




1.2 X 


10 


-4^1/2^-1/2^-1 


2.1 


X 


10" 




Ion electron 


















■tie/ ry2 ^. 

I/g' /UeZ Xie 




1.6 X 


10 


-9^-1^-3/2 


1.7 


X 


10" 


4 1/2 -3/2 

jji ' e ' 


/UeZ Xie 




3.2 X 


10 


-9^-1^-1/2^-1 


1.8 


X 


10" 


7^-1/2^-3/2 


i\e / ryl \ 

V ' jUeZ Xie 




1.6 X 


10 


-9^-1^-1/2^-1 


1.7 


X 


10" 


V'/'Te"^/' 



Ion— ion 



/1/2 / / \ -1/2 

-6.8X10"^^^^ 1 T-^/^ 



-> 9.0 X 10 - + — ^-TTT 



i\i' 



n^,Z^Z'^X^^, 



1.4xlO-V'/V"'T-^/^e-^ 

^ 1.8 X 10" V'/'e"^/^ 

6.8 X 10" V'/'^-'T-^/^e"^ 

^9.0 X lO-V'^V'-'Te"^/' 
In the same limits, the energy transfer rate follows from the identity 

i/g = 2vs - u± - u\\ , 

except for the case of fast electrons or fast ions scattered by ions, where the 
leading terms cancel. Then the appropriate forms are 

uf^ — y 4.2 X lO'^mZ'^Xei 



"^/^/x"^ -8.9 X 10^(Ai/T)^/^e"^ exp(-1836Aie/T) 



-1 

sec 



32 



and 

1^*1'' — y 1.8 X 10~^ n^^Z^Z'^X^^, 

In general, the energy transfer rate positive for e > Cq,* and nega- 

tive for e < €0,*, where cc* = {mp /mct)£a* /Tp is the solution of ip' (x*) = 
(mo; |m;3)'?/;(cc*). The ratio €01* /Tp is given for a number of specific a, /3 in the 
following table: 



a\/3 




e e 


e, i 


i e\p e 


D e 


T, e 


He^ e 


He^ 


f * 


1.5 


0.98 


4.8 X 10"^ 2.6 X 


10"^ 1.8 X 10"^ 1.4 X 10"^ 



When both species are near Maxwellian, with Tj ^ Tg, there are just 
two characteristic collision rates. For Z — 1, 

Ue = 2.9 X 10"®nATe"^/^ sec"^ 

Ui =4.8 X 10"®nATi"^/^//"^/^ sec"\ 



Temperature Isotropization 

Isotropization is described by 



dt 



1 dTn 



2 dt 



where, if A = T^/Tii - 1 > 0, 



2 2 

cx^ap A — 2 



-3 + (A + 3) 



tan-^(A^/") 



If A < 0, tan-^(A^/2)/A^/2 is replaced by tanh-^(-A)^/2/(-A)^/2. For 
T|| = T, 



i/^ = 8.2 X 10"'^nAT"^/^ sec"^ 

i/^ = 1.9 X 10"®nAZ^//"^/^T"^/^ sec 
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Thermal Equilibration 

If the components of a plasma have different temperatures, but no rela- 
tive drift, equilibration is described by 



dt 



where 

^ = 1.8 X 10 — sec 

{maT/3 + mj3To,Y'^ 

For electrons and ions with K. Ti = T, this implies 

uf'/m = ui^yue = 3.2 X 10~^Z^A///T^/^cm^ sec~\ 

Coulomb Logarithm 

For test particles of mass rria. and charge Cq, = Zae scattering off field 
particles of mass and charge = Zf^e, the Coulomb logarithm is defined 

as A = In A = ln(rmax/^min)- Here r„iin is the larger of eQ,e^/mQ,^u^ and 
h/2mctf3U, averaged over both particle velocity distributions, where mcef3 = 
mcm/3 /{rric + m/3) and u = Vq, - V/3; r^ax = (47r ^ n^e^^/fcT^)"-*-/^, where 
the summation extends over all species 7 for which < VT-y'^ = kT^/rrij. If 

this inequality cannot be satisfied, or if either uu!ca~^ < ?"max or uuJci3~^ < 
fmaic, the theory breaks down. Typically A » 10—20. Corrections to the trans- 
port coefficients are 0(A~^); hence the theory is good only to ~ 10% and fails 

when A ~ 1. 

The following cases are of particular interest: 

(a) Thermal electron— electron collisions 

Aee = 23.5 - ln(ne'/'Te-"/^) - [10-^ + (InTe - 2)Vl6]'/' 

(b) Electron— ion collisions 



Arx — A^ 



23 



— In (jie^ 



/2 



ZT' 



-3/2 



Tim^/mi < Te < lOZ eV; 
= 24-ln (ne^/^T"^), Time/rm < lOZ^ eV < Te 

= 30 - In (ni^^^Ti~^^^Z^fj,~^y Tg < TiZme/rrii. 

(c) Mixed ion— ion collisions 



^ii' — ^i'i — 23 



In 



z z' {jji + jjb') f mZ^ ^ '^iiZ 



12 \ 1/2 



+ 
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(d) Counterstreaming ions (relative velocity vd = Pdc) in the presence of 
warm electrons, kTi/rrii, kT^/ /m^/ < vn^ < kTe/rUe 



\i' = \'i = 35 - In 



Fokker- Planck Equation 



Dt 



dt 



+ V V/" + F- Vv/' 



or 

dt 



coll 



where F is an external force field. The general form of the collision integral is 



run, 



d%' (u^ I — uu)u ^ 

|— r(v)v^,/^(v') - J-/^(v')Vvr(v)| 



(Landau form) where u = v' — v and / is the unit dyad, or alternatively. 



J"\^ = 47rA«/3^^^ |r (v)Vvi/(v) - ■ [r (v)VvVvG(v)] | , 



rrif 



where the Rosenbluth potentials are 



G(v) = / f^{-v')ud\' 



H(v) = 1 + 



m/3 




nl3 t l\ -1 ,3 / 

J (y )u a V . 



If species a is a weak beam (number and energy density small compared with 
background) streaming through a Maxwellian plasma, then 



tQ;\/3 "^a a\/3 j-a. ^ a\/3 ^7 
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B-G-K Collision Operator 

For distribution functions with no large gradients in velocity space, the 
Fokker-Planck collision terms can be approximated according to 



Dt 



= J^ee{Fe - fe) + l^ei{Fe - fe)\ 



Dfi 
Dt 



= J^ie{Fi - fi) + Vii{Fi - fi). 



The respective slowing-down rates v^^^ given in the Relaxation Rate section 
above can be used for i^a/g, assuming slow ions and fast electrons, with e re- 
placed by Tq,. (For u^e and ua, one can equally well use and the result 
is insensitive to whether the slow- or fast-test-particle limit is employed.) The 
Maxwellians Fa and Fa are given by 



'a — 



(■ 



rric 



3/2 



exp 



ma{-v - Vg) 



2 n 



Fa — Tla ' 



m, 



3/2 



V 27rkTa 



exp < - 



n2 



2kTn 



where 

5 ^ a and Ta are the number density, mean drift velocity, and effective 
temperature obtained by taking moments of . Some latitude in the definition 
of Ta and Vq, is possible;^° one choice is Tg = Ti, Ti = Te, Vg = v^, = Vg. 

Transport Coefficients 

Transport equations for a multispecies plasma: 



d'^ria 
dt 



+ no, V • Vq, = 0; 



dt 



= -Vpa - V ■ Pa+Za 



en. 



E + -v« X B 

c 



3 d k 'J. 

— "Tla 

2 



dt 



+ Pa V • Vq, = -V • Qq - Pq : VVq, + Qc 



Here d"/dt = d/dt + 

■ V; Pq, — HakTai where k is Boltzmann's constant; 

Rq, = and Qa = Qa/3, where Rq,/3 and Qa/s are respectively 

the momentum and energy gained by the ath species through collisions with 
the /3th; Pa is the stress tensor; and Qq, is the heat flow. 
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The transport coefficients in a simple two-component plasma (electrons 
and singly charged ions) are tabulated below. Here || and _L refer to the di- 
rection of the magnetic field B = hB; u = Ve — is the relative streaming 
velocity; Ue = rii = n; j = — neu is the current; ujce = 1-76 X lO^Ssec""*^ and 
<-^ci = {iTie/ iTii)u}ce are the electron and ion gyrofrequencies, respectively; and 
the basic collisional times are taken to be 



^ ^ ^ = 3.44 X 10*^ 



rjp 3/2 



4V27rnAe4 



sec, 



where A is the Coulomb logarithm, and 



Ti = 



4V7rn Ae'^ 



3/2 



2.09 X 10 



3/2 

7 -'-i 1/2 



nX 



11 sec. 



In the limit of large fields {ujcaTot a. = i,e) the transport processes may 

be summarized as follows: 



momentum transfer 
frictional force 

electrical 
conductivities 

thermal force 

ion heating 

electron heating 
ion heat flux 

ion thermal 
conductivities 
electron heat flux 

frictional heat flux 

thermal gradient 
heat flux 

electron thermal 
conductivities 



R-ei — "R-ie = R — R-u + R-T ; 
Ru = ne(j||/(T|| +j±/cT±); 

cr|| = 1.96(T_l; a± = ne Te/irie', 

3n 

R^= -0.71nV|| (fcTe) b X V ±(kTe); 



SrUp nk 



— (Te-Ti); 



Qi 

Qe 

c[i = -^V||(fcTi) - k'^V ±(kTi) + K\h X V^(fcTi); 



rrii Te 
-Qi — R • u; 



Kn — 3.9- 



nkTiTi 
rrii 

qe ^Qu + Qt; 



2nkTi 
rriiujjri ' 



— 



5nkTi 



e 



3nkT 

0.71nfcTeU|| H x u^; 

-KjV||(fcTe) - KlV^(fcTe) - K^h X V^(fcTe); 



= 3.2- 



nkTpTp 



4.7- 



nkTf. 



rrif 



2, 
' ce 



A 



5nkTe 

2me(jJc 
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stress tensor (either Pxx= — — {W^x + ^yy) —{Wxx 

species) 2 2 



^yy) - m^xy\ 



Pyy= -"^{Wxx+Wyy) + !^l(Wxx-Wyy) + mWxy; 

773 

Pxy= Pyx = -rjl^xy + —(Wxx " Wyy)] 



Pxz = Pzx = —m^xz - ri4^Wyz] 
Pyz= Pzy = -m^yz + 'nAWxz] 

Pzz = -rjoWzz 
(here the z axis is defined parallel to B); 

SnkTi 



ion viscosity 



electron viscosity 



I 

V3 



e 
^0 



e 
^3 



= 0.96nkTiTi; r){ = 



10a; ^Ti 



nkT, 



2oj, 



I i 

-; V4 = 



nkTi 



= 0.73nkTeTe; ill - 0.51 



n/cTe 

2Ur 



; ^4 = 



nkTe 



'ce ^ce 

For both species the rate- of- strain tensor is defined as 



nkT^ 



QnkTi 
5a; ?ri ' 



; r7| = 2.0^^, 



w. 



dv-j dv 



jk 



+ 



k 



dxk dxj 3 



- -Sjk^ ■ V. 



When B = the following simplifications occur: 



Ru = nej/c7||; Rt = -0.71nV(fcTe); = -k?, V(fcTi); 



= O.rinfcTeU; q^ = -KjV(fcTe); Pjfe = -T7oWjfe- 

For ojceTe 1 ^ ujciTi, the electrons obey the high-field expressions and the 
ions obey the zero-field expressions. 

Collisional transport theory is applicable when (1) macroscopic time rates 
of change satisfy d/dt <C where r is the longest collisional time scale, and 

(in the absence of a magnetic field) (2) macroscopic length scales L satisfy L ^ 
I, where I = vt is the mean free path. In a strong field, ujce^ ^ I5 condition 
(2) is replaced by ^ I and L± ^ VlrZ {L± ^ re in a uniform field), 
where Ly is a macroscopic scale parallel to the field B and L± is the smaller 
of -B/|V_L-B| and the transverse plasma dimension. In addition, the standard 
transport coefficients are valid only when (3) the Coulomb logarithm satisfies 

A ^ 1; (4) the electron gyroradius satisfies re ^ Ad, or STrneTrieC^ ^ B^; (5) 
relative drifts u = Vq, — V/3 between two species are small compared with the 
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thermal velocities, i.e., <^ kTa/ma, kTjj/mj^; and (6) anomalous transport 
processes owing to microinstabilities are negligible. 

Weakly Ionized Plasmas 

Collision frequency for scattering of charged particles of species a. by 
neutrals is 

a|0/; rn / \l/2 

where no is the neutral density, cr"^*^ is the cross section, typically 5 X lO""*"^ 

cm^ and weakly dependent on temperature, and (To/mo)''^^^ < (Ta/ma)^^'^ 
where Tq and mo are the temperature and mass of the neutrals. 

When the system is small compared with a Debye length, L <C ^d, the 
charged particle diffusion coefficients are 

Dcx. — kT (X I m oi^v ^ 

In the opposite limit, both species diffuse at the ambipolar rate 



Da = 



where /Xq; = ea/mcti^cx is the mobility. The conductivity a a satisfies ctq, = 
In the presence of a magnetic field B the scalars /x and cr become tensors, 

J" = • E = aJ^Eii + (T^E^ + (T^E X b, 
where h — B and 

a ot I ( 2. 2\ 

OL oi / / 2 I 2 \ 

CT/^ = cry i/ai^ca/ii^a + ^ coc) ■ 

Here a±_ and cr^ are the Pedersen and Hall conductivities, respectively. 
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APPROXIMATE MAGNITUDES 
IN SOME TYPICAL PLASMAS 



Plasma Type 


3 

n cm 


1 e V 




sec 


A D cm 




sec ^ 


Interstellar gas 


1 


1 


6 


X 


10^ 


7 


X 


10^ 




4 


X 10^ 


7 


X 10"^ 


Gaseous nebula 


10^ 


1 


2 


X 


10^ 




20 




8 


X 10^ 


6 


X 10"^ 


Solar Corona 


10^ 


10^ 


2 


X 


10^ 


2 


X 


10" 


1 


8 


X 10^ 




60 


uiiiusti not piabnia 


lU 


1 

lU 


6 


X 


lU 


7 


X 


10" 


3 


4 


X 10^ 




40 


Solar atmosphere, 


lU 


1 


6 


X 


mil 


7 


X 


10" 


5 




40 


2 


X 10^ 


gas discharge 




























Warm plasma 




10 


6 


X 




2 


X 


10" 


4 


8 


X 10^ 




10^ 


Hot plasma 




10^ 


6 


X 




7 


X 


10" 


4 


4 


X 10^ 


4 


X 10^ 


Thermonuclear 




10^ 


2 


X 




2 


X 


10" 


3 


8 


X 10^ 


5 


X 10^ 


plasma 




























Theta pinch 




10^ 


6 


X 




7 


X 


10" 


5 


4 


X 10^ 


3 


X 10^ 


Dense hot plasma 




10^ 


6 


X 




7 


X 


10" 


6 


4 


X 10^ 


2 


X 10^0 


Laser Plasma 


10^0 


10^ 


6 


X 


10^4 


7 


X 


10" 


7 




40 


2 


X 10^^ 



The diagram (facing) gives comparable information in graphical form. 



40 



25 



50% ionization 
hydrogen plasma 



20 - 



15 - 



CO 



6 

o 



O 



10 



(4TT/3)nXjj < 1 



Laser 
Plasma 



Focus 





Shock 
tubes 






Z— pinches 


/ High 










/ pressure 






/ arcs 







Alkali 
metal, '' 
plgisma 



Flames 



5 - 



Earth 
iono- 
sphere 







Low 
pressure 



Fusion 
experiments 



Glow 
discharge 



Solar 
corona 



> 1 cm 



Solar wind 
(1 AU) 



Earth 
plasma sheet 







1 2 
log^O T (ev) 



A. 

1- 

A 



> 1 cm 



Fusion 
reactor 



4 



5 
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IONOSPHERIC PARAMETERS 



The following tables give average nighttime values. Where two numbers 
are entered, the first refers to the lower and the second to the upper portion 
of the layer. 



Quantity 


E Region 


F Region 


AltitiiHp /^km'l 


Qf)— 1 fin 

kJ\J -L \J\J 


1 fif)— 500 


Number density (m ) 


1.5 X 10 -3.0 X 10 


5 X 10 -2 X 10 


Height-integrated number 


9 X lO^'* 


4.5 X 10 


density (m~ ) 






T J 1 11"" 

Ion-neutral collision 


2 X 10-10 


0.5-0.05 


frequency (sec ) 






T / n • • 

Ion gyro- /collision 


0.09-2.0 


4.6 X 10 -5.0 X 10 


frequency ratio Ki 






Ion Pederson factor 


0.09-0.5 


2.2 X 10 -2 X 10 


/ / 1 1 2 \ 
Ki/{1 + Ki ) 






T T T 1 1 J? J_ 

Ion Hall ractor 


8 X 10 -0.8 


1.0 


i^x / rit ) 






. 1 if*{~\T\ — nciii'i'T'Ql /^olliGion 


1 ^ V 1 (T*— Q n V 1 n-^ 


80—1 n 


f r pn 1 1 pn cv 

1.1. vI^VJ LAt/XX^ y 






Electron gyro-/collision 


4.1 X 10^-6.9 X 10^ 


7.8 X 10^-6.2 X 10^ 


frequency ratio Ke 






Electron Pedersen factor 


2.7 X 10"^-1.5 X 10"^ 


10"^-1.5 X 10"^ 








Electron Hall factor 


1.0 


1.0 








Mean molecular weight 


28-26 


22 16 


Ion gyrofrequency (sec""*^) 


180-190 


230-300 


Neutral diffusion 


30-5 X 10^ 


10^ 


coefficient (m^sec""*^) 







The terrestrial magnetic field in the lower ionosphere at equatorial latti- 
tudes is approximately Bq = 0.35 x lO"'* tesla. The earth's radius is Re = 6371 
km. 
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SOLAR PHYSICS PARAMETERS 



Parameter 


Symbol 


Value 


Units 


Total mass 


Mq 


1.99 X 10^^ 


g 


Radius 


Rq 


6.96 X 10^° 


cm 


Surface gravity 


9Q 


2.74 X 10^ 


cm s 


Escape speed 




6.18 X 10^ 


cm s 


Upward mass flux in spicules 





1.6 X 10"^ 


2 1 

gem s 


Vertically integrated atmospheric density 





4.28 


gcm~^ 


Sunspot magnetic fleld strength 


-Smax 


2500-3500 


G 


Surface effective temperature 


To 


5770 


K 


Radiant power 




3.83 X 10^^ 


ergs 


Radiant flux density 


JF 


6.28 X 10^° 


ergcm~^s~^ 


Optical depth at 500 nm, measured 




0.99 




from photosphere 








Astronomical unit (radius of earth's orbit) 


AU 


1.50 X 10^^ 


cm 


Solar constant (intensity at 1 AU) 


/ 


1.36 X 10^ 


erg cm~^ s~ 



Chromosphere and Corona 



Parameter (Units) 


Quiet 
Sun 


Coronal 
Hole 


Active 
Region 


Chromospheric radiation losses 










(ergcm~^ s~^) 










Low chromosphere 




2 X 10*^ 


2 X 10^ 


> 10^ 


Middle chromosphere 




2 X 10*^ 


2 X 10*^ 


10" 


Upper chromosphere 




3 X 10^ 


3 X 10^ 


2 X 10^ 


Total 




4 X 10*^ 


4 X 10*^ 


> 2 X lO" 


Transition layer pressure (dyne cm 




0.2 


0.07 


2 


Coronal temperature (K) at 1.1 Rq 


1.1-1.6 X 10^ 


10^ 


2.5 X 10^ 


Coronal energy losses (ergcm^'^s" 










Conduction 




2 X 10^ 


6 X 10^ 


lO^-lO" 


Radiation 




10^5 


10^ 


5 X 10^ 


Solar Wind 




< 5 X 10^ 


7 X 10^ 


< 10^ 


Total 




3 X 10^ 


8 X 10^ 


10" 


Solar wind mass loss (gcm"'^ s""*^) 




< 2 X 10"^^ 


2 X 10"^° 


< 4 X 10"^^ 
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THERMONUCLEAR FUSION^^ 

Natural abundance of isotopes: 

hydrogen '^nl'^H — ^-^ ^ 10~^ 

helium ^He^/^He^ ~ 1-3 X 10~® 

lithium nLie/^T-Li? = 0.08 

Mass ratios: rrie/mD = 2.72 x 10~^ = 1/3670 

(me/mD)^/^= 1.65 x 10"^ = 1/60.6 
me/rriT = 1.82 X 10"^ = 1/5496 

(me/mT)^/^= 1.35 X 10"^ - 1/74.1 

Absorbed radiation dose is measured in rads: 1 rad = lO'^ erg g""*^ . The curie 
(abbreviated Ci) is a measure of radioactivity: 1 curie = 3.7 X 10^° counts sec~^ . 

Fusion reactions (branching ratios are correct for energies near the cross section 
peaks; a negative yield means the reaction is endothermic): 

^T(l.OlMeV) + p(3.02MeV) 
>He^(0.82MeV) + n(2.45MeV) 
>He^(3.5MeV) + n(14.1MeV) 
>He^(3.6MeV) + p(14.7MeV) 



(la) 


D + D - 


50% 


(lb) 




50% 


(2) 


D + T - 




(3) 


D + He^- 




(4) 


T + T - 




(5a) 


He^ + T- 


51% 


(5b) 




43% 


(5c) 


P + Li^ - 


6% 


(6) 




(7a) 


P + Li^ - 


20% 


(7b) 


D + Li^ - 


80% 


(8) 




(9) 


P + B^^- 




(10) 


n + Li^ - 





>He^ + 2n + 11.3 MeV 



>He^ + p + n + 12.1 MeV 
>He^(4.8MeV) + D(9.5 MeV) 



>He^(1.7MeV) + He^(2.3MeV) 



>2 He^ + 17.3 MeV 



^Be"^ + n - 1.6 MeV 



>2He^ + 22.4 MeV 



>3 He^ + 8.7 MeV 



>He^(2.1MeV) + T(2.7MeV) 



The total cross section in barns (1 barn = 10~^* cm^) as a function of E, the 
energy in keV of the incident particle [the first ion on the left side of Eqs. 
(1)— (5)], assuming the target ion at rest, can be fitted hy^^°' 

A5+ \(A4-A3Ef + iyA2 
ctt{E) = - 



E [exp(Ai£;-i/2) _ ij 
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where the Duane coefficients Aj for the principal fusion reactions are as follows: 





D D 


D D 


D-T 


D-He^ 


T T 


T He"^ 




(la) 


(lb) 


(2) 


(3) 


(4) 


(5a-c) 




46.097 


47.88 


45.95 


89.27 


38.39 


123.1 


A2 


372 


482 


50200 


25900 


448 


11250 


A3 


4.36 X 10-"^ 


3.08 X 10"* 


1.368 X 10-2 


3.98 X 10"^ 


1.02 X 10-3 







1.220 


1.177 


1.076 


1.297 


2.09 





As 








409 


647 









Reaction rates av (in cm^ sec averaged over Maxwellian distributions: 



Temperature 




D 


-D 






D 


-T 




D-He^ 






T 


-T 




T-He^ 




(kcV) 


(la- 


f lb) 




(2) 




(3) 






(4) 




(5a-c) 




1.0 


1.5 


X 


10" 


22 


5.5 


X 


10" 


21 


lQ-26 




3.3 


X 


10" 


22 


lQ-28 




2.0 


5.4 


X 


10" 


21 


2.6 


X 


10" 


19 


1.4 X 10" 


23 


7.1 


X 


10" 


21 


lQ-25 




5.0 


1.8 


X 


10" 


19 


1.3 


X 


10" 


17 


6.7 X 10" 


21 


1.4 


X 


10" 


19 


2.1 X 10" 


22 


10.0 


1.2 


X 


10" 


18 


1.1 


X 


10" 


16 


2.3 X 10" 


19 


7.2 


X 


10" 


19 


1.2 X 10" 


20 


20.0 


5.2 


X 


10" 


18 


4.2 


X 


10" 


16 


3.8 X 10" 


18 


2.5 


X 


10" 


18 


2.6 X 10" 


19 


50.0 


2.1 


X 


10" 


17 


8.7 


X 


10" 


16 


5.4 X 10" 


17 


8.7 


X 


10" 


18 


5.3 X 10" 


18 


100.0 


4.5 


X 


10" 


17 


8.5 


X 


10" 


16 


1.6 X 10" 


16 


1.9 


X 


10" 


17 


2.7 X 10" 


17 


200.0 


8.8 


X 


10" 


17 


6.3 


X 


10" 


16 


2.4 X 10" 


16 


4.2 


X 


10" 


17 


9.2 X 10" 


17 


500.0 


1.8 


X 


10" 


16 


3.7 


X 


10" 


16 


2.3 X 10" 


16 


8.4 


X 


10" 


17 


2.9 X 10" 


16 


1000.0 


2.2 


X 


10" 


16 


2.7 


X 


10" 


16 


1.8 X 10" 


16 


8.0 


X 


10" 


17 


5.2 X 10" 


16 



For low energies (T ^ 25 keV) the data may be represented by 

(o^)oD =2.33 X 10"^^T"2/3exp(-18.76T"^/3)cm3sec"^ 



(av)z)T = 3.68 X 10 ^^T exp(-19.94T ^/^^ ^,^^^3 1^ 
where T is measured in keV. 

A three-parameter model has also been developed for fusion cross-sections of 
hght nuclei. 

The power density released in the form of charged particles is 

Pdd = 3.3 X 10~^^nn^('av)DDWdittcni~^ (including the subsequent 

prompt D— T reaction only) ; 

Pdt = 5.6 X 10~^^nD'nTioU) DT 'Wdittcm~^; 

^DHe3 — 2-9 X 10~^^'^Z3'^He3 (^)r»He3 wattcm"^. 
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RELATIVISTIC ELECTRON BEAMS 

Here 7 = (1 — is the relativistic scaling factor; quantities in 

analytic formulas are expressed in SI or cgs units, as indicated; in numerical 
formulas, / is in amperes (A), B is in gauss (G), electron linear density N is 
in cm~^, and temperature, voltage and energy are in MeV; jSz = Vz/c; k is 
Boltzmann's constant. 

Relativistic electron gyroradius: 

(7^ - 1)^/^ (cgs) = 1.70 X 10^(7^ - 1)^/^5"^ cm. 



fe = 

eB 



Relativistic electron energy: 



W - rac-i = O.5II7 MeV. 



Bennett pinch condition: 



= 2Nk{Te + T^c^ (cgs) = 3.20 X 10"^Ar(Te + A^. 
Alfven-Lawson limit: 

Ia = (mcVe)^z7 (cgs) = (47rmc//xoe)/3^7 (SI) = 1.70 x 10^/3^7 A. 
The ratio of net current to I a is 

/ _ 1/ 
Ia 7 

Here u = Nr^ is the Budker number, where re = e /mc = 2.82 x 10 cm 
is the classical electron radius. Beam electron number density is 

Ub = 2.08 X 10^JP~^ cm"^, 

where J is the current density in Acm~^. For a uniform beam of radius a (in 
cm), 

Ub = 6.63 X 10^/a~^/3~^ cm~^, 



and 



2re 



a 7 
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Child's law: (non-relativistic) space-charge-limited current density between 
parallel plates with voltage drop V (in MV) and separation d (in cm) is 



J = 2.34 X lO^y^/^d"^ Acm"^. 



The saturated parapotential current (magnetically self-limited flow along equi- 
potentials in pinched diodes and transmission lines) is^^ 

Ip = 8.5 X 10^G7 In 7 + (7^ - 1)^^^] A, 

where G is a geometrical factor depending on the diode structure: 

^ w for parallel plane cathode and anode 

27rd of width w, separation d; 

G = ^In cylinders of radii R\ (inner) and R2 (outer); 

^ Rc for conical cathode of radius Rc-, maximum 

do separation do (at r = Rc) from plane anode. 

For /3 — >■ (7 — >■ 1), both I a and Ip vanish. 

The condition for a longitudinal magnetic fleld to suppress fllamentation 
in a beam of current density J (in Acm~^) is 

> 47/3, (7 J)'/^ G. 

Voltage registered by Rogowski coil of minor cross-sectional area A, n turns, 
major radius a, inductance L, external resistance R and capacitance C (all in 
SI): 

externally integrated V = {1/ RC){nAiXQl /2tt a); 

self-integrating V = {R/ L)(nAiJ,oI /27Ta) = Rl/n. 

X-ray production, target with average atomic number Z (y < 5 MeV): 
?7 = x-ray power/beam power = 7 X 10 ^ ZV. 

X-ray dose at 1 meter generated by an e-beam depositing total charge Q 
coulombs while V > 0.84V'inax in material with charge state Z: 

D = IdOVjjQZ'^^ rads. 
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BEAM INSTABILITIES 



Name 


Conditions 


Saturation Mechanism 


Electron- 
electron 


Vd > Ve,-, j = 1,2 


Electron trapping until 

Ve,- ~ Vd 


Buneman 


Vd > {M/my/^Vi, 

Vd > Ve 


Electron trapping until 

Ve ~ Vd 


Beam-plasma 


Vb > (np/nf,)'/^V6 


Trapping of beam electrons 


Weak beam- 
plasma 


Vb < inp/nty/^Vb 


Quasilinear or nonlinear 
(mode coupling) 


Beam-plasma 
(hot-electron) 


Ve > n > n 


Quasilinear or nonlinear 


Ion acoustic 


Te > Ti, Vd > Cs 


Quasilinear, ion tail form- 
ation, nonlinear scattering, 
or resonance broadening. 


Anisotropic 

lt3Ill]Jt31 clL HI c 

(hydro) 


Te_L > 2Te|| 


Isotropization 


Ion cyclotron 


T/, \ oaf/ (t^-r- 

Vd > ^yjVi \Jor 


Ion heating 


JDcclIIl-CyClOljI oil 

(hydro) 


Vd > 


rtcooiiancc uioaQciiiiig 


ivioainea xwo- 
stream (hydro) 


y d <■ \i- -r P ) V Ai 
Vd > Cs 


Trapping 


Ion-ion (equal 
beams) 


U < 2(1 + /3)1/Va 


Ion trapping 


Ion-ion (equal 
beams) 


U < 2Cs 


Ion trapping 



For nomenclature, see p. 50. 
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Name 



Parameters of Most Unstable Mode 



Growth Rate 



Frequency 



Wave 
Number 



Group 
Velocity 



Electron- 
electron 

Buneman 
Beam-plasma 



Weak beam- 
plasma 

Beam-plasma 
(hot-electron) 

Ion acoustic 

Anisotropic 
temperature 
(hydro) 

Ion cyclotron 

Beam-cyclotron 
(hydro) 

Modified two- 
stream 
(hydro) 

Ion-ion (equal 
beams) 

Ion-ion (equal 
beams) 



— OJe 

2 



/ m 
0.7 — 

0.7 



1/3 



1/3 



CO, 



rib ( H 



2np \Vb 



m 



M 



1/2 



(jOi 



He 

O.lQi 
0.7Qe 



1 



OAuJi 







0.4 — uj. 



00 e — 

0.4 — UOe 
\npj 



00, 



^—00^ 
Ve 



OOi 

OOe COS 6 Q,^ 

1.2Qi 






OOe 

0.9 — 

Vd 



OJe 
OJe 



OJe 



1.7- 



H 



1.2 



Vd 

u 



OJi 

1.2 — 

U 







2 
2 

3^^ 



Ve 



1 - 

> Vd; 
< 
1 








For nomenclature, see p. 50. 
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In the preceding tables, subscripts e, i, d, b, p stand for "electron," "ion," 
"drift," "beam," and "plasma," respectively. Thermal velocities are denoted 
by a bar. In addition, the following are used: 



m 


electron mass 




gyroradius 


M 


ion mass 




plasma/magnetic energy 


V 


velocity 




density ratio 


T 


temperature 


Va 


Alfven speed 




number density 




gyrofrequency 


n 


harmonic number 




hybrid gyrofrequency. 


Cs = (Te/M)^/2 


ion sound speed 








plasma frequency 


U 


relative drift velocity of 


Xd 


Debye length 




two ion species 
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LASERS 



System Parameters 

Efficiencies and power levels are approximate.^ 



Type 


Wavelength 


Hjinciency 


Power levels available (W) 


(^xm) 


Jr uisccl 


L- VV 


CO2 


10.6 


0.01-0.02 
(pulsed) 


> 2 X 10^3 


> 10^ 







U.4 


> 10 


> 100 


Holmium 


z.Uo 


U.UoJ— U.iI|I 


> 10 


80 


ioame 


i.oio 


U.UUo 


3 X 10 




IN a-giass 


1 .UD 




1.25 X 10 




1\T A . "V A 

IN a: Y AVor 


i.UD4 




-1 n9 
10 


> 10 


In Q. I Lir 


1 n/i 

i .U41:0, 

1.54,1.313 




4 X 10 


80 


Nd:YV04 


1.064 






> 20 


Er:YAG 


2.94 


_ 


1.5 X 10^ 




*Color center 


1-4 


10~^ 


5 X 10^ 


1 


*Ti:Sapphire 


0.7-1.5 


0.4 X r]p 


10^4 


150 


Ruby 


0.6943 


< 10-3 


10^0 


1 


He-Ne 


0.6328 


10-^ 


— 


1-50x10-3 


*Argon ion 


0.45-0.60 


10-3 


5 X 10* 


150 


*OPO 


0.3-10 


> 0.1 X rjp 


10^0 


5 


N2 


0.3371 


0.001-0.05 


10® 




*Dye 


0.3-1.1 


10-3 


5 X 10^ 


> 100 


Kr-F 


0.26 


0.08 


10^2 


500 


Xenon 


0.175 


0.02 


> 10^ 




Ytterbium fiber 


1.05-1.1 


0.55 


5 X 10"^ 


10* 


Erbium fiber 


1.534 




7 X 10® 


100 


Semiconductor 


0.375-1.9 


> 0.5 


3 X 10^ 


> 103 


*Tunable sources 


f lamp-driven 


^diode-driven 





Nd stands for Neodymium; Er stands for Erbium; Ti stands for Titanium; 
YAG stands for Yttrium— Aluminum Garnet; YLF stands for Yttrium Lithium 
Fluoride; YV05 stands for Yttrium Vanadate; OPO for Optical Parametric 
Oscillator; r]p is pump laser efficiency. 
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Formulas 

An e-m wave with k || B has an index of refraction given by 

where ± refers to the hehcity. The rate of change of polarization angle as a 
function of displacement s (Faraday rotation) is given by 

d9/ds = {k/2){n- - n+) = 2.36 X lO^ATB/"^ cm"\ 

where N is the electron number density, B is the field strength, and / is the 
wave frequency, all in cgs. 

The quiver velocity of an electron in an e-m field of angular frequency uj 

is 

Vq = eEma.x/'mi^ — 25.6/ ' Aq cm sec 

in terms of the laser fiux I = cE^^^/Stt, with I in watt/cm^, laser wavelength 
Ao in /xm. The ratio of quiver energy to thermal energy is 

Wqu/Wth = m^vo^/2kT = 1.81 x IQ-^^Ao^i/T, 

where T is given in eV. For example, if / = lO^^Wcm"'^, Aq = 1 yiim, T = 
2keV, then W^^/Wt\, ^0.1. 

Pondermotive force: 

= A/'V(£;^)/87rA/"c, 

where 

= 1.1 X 10^Uo~^cm"^. 

For uniform illumination of a lens with /-number F, the diameter d at 
focus (85% of the energy) and the depth of focus / (distance to first zero in 
intensity) are given by 

d K 2AAF\e/9DL and I ^ ±2F'^\e/enL- 

Here 6 is the beam divergence containing 85% of energy and Odl is the 
diffraction-limited divergence: 

9dl = 2.44A/6, 

where b is the aperture. These formulas are modified for nonuniform (such as 
Gaussian) illumination of the lens or for pathological laser profiles. 
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ATOMIC PHYSICS AND RADIATION 



Energies and temperatures are in eV; all other units are cgs except where 
noted. Z is the charge state (Z = refers to a neutral atom); the subscript e 
labels electrons. N refers to number density, n to principal quantum number. 
Asterisk superscripts on level population densities denote local thermodynamic 
equilibrium (LTE) values. Thus AT^* is the LTE number density of atoms (or 
ions) in level n. 

Characteristic atomic collision cross section: 

(1) Trao^ = 8.80 X 10"^'^cm^. 

Binding energy of outer electron in level labelled by quantum numbers n, Z: 



(2) E^(n,l) = -- 

where = 13.6 eV is the hydrogen ionization energy and — 0.75Z , 

I ^ 5, is the quantum defect. 

Excitation and Decay 

Cross section (Bethe approximation) for electron excitation by dipole 
allowed transition m — >■ n (Refs. 32, 33): 

(3) amn = 2.36 X iq-13 ^) ,^2^ 

where fmn is the oscillator strength, g(n,m) is the Gaunt factor, e is the 
incident electron energy, and AEnm — En — Em- 

Electron excitation rate averaged over Maxwellian velocity distribution, X^n 
= Neiamnv) (Refs. 34, 35): 



(4) Xmn = 1-6 X 10 — — exp -- 



nm 



'-^^nm ^ e 



where {g{n, m)) denotes the thermal averaged Gaunt factor (generally ~ 1 for 
atoms, ^ 0.2 for ions). 
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Rate for electron collisional deexcitation: 



(5) y^nm — {Njn* / Nn*)Xmn- 

Here Nm* / Nn* = (gm / dn) ^^pi^Enm /Te) is the Boltzmann relation for level 
population densities, where gn is the statistical weight of level n. 

Rate for spontaneous decay n ^ m (Einstein A coefficient)^* 

(6) Anm = 4.3 X 10 {gm / gn) fmni^Enm) SeC . 



Intensity emitted per unit volume from the transition n — >■ m in an optically 
thin plasma: 

(7) Inm = 1.6 X 10" ^^A^^AT^A^^^ watt/cm^. 



Condition for steady state in a corona model: 

(8) NoN^{aor,v) ^ Nr^Ar^O, 

where the ground state is labelled by a zero subscript. 
Hence for a transition n — >■ m in ions, where {g{n, 0)) K. 0.2, 

3 



-25 fnmgmNeNo /AEnm\ f AEno\ Watt 



cm"^ 



Ionization and Recombination 

In a general time-dependent situation the number density of the charge 
state Z satisfies 



. , dN(Z) 

(10) -ir = ^' 



- S{Z)N(Z) - a{Z)N{Z) 
+S(Z - 1)N(Z - 1) + a(Z + 1)N(Z + 1) 



Here S{oZ) is the ionization rate. The recombination rate (y:{Z) has the form 
a(Z) = oir{Z) + Nf>OLs{Z), where otr and 0:3 are the radiative and three-body 
recombination rates, respectively. 
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Classical ionization cross-section^^ for any atomic shell j 

(11) C7i = 6 X 10~^'^bjgj{x)/Uj'^ cni^. 

Here bj is the number of shell electrons; Uj is the binding energy of the ejected 

electron; x = e/Uj, where e is the incident electron energy; and is a universal 
function with a minimum value g^min ~ 0.2 at cc ~ 4. 

Ionization from ion ground state, averaged over Maxwellian electron distribu- 
tion, for 0.02 < Te/E^ < 100 (Ref. 35): 

(12) S(Z) = 10"^ . ^.oSo °° exp - ^ cmVsec, 

where E^ is the ionization energy. 

Electron-ion radiative recombination rate (e -|- N(Z) — >■ N(Z — 1) -|- hv) 
for Te/Z^ < 400 eV (Ref. 37): 



(13) ar(Z) = 5.2 X 10"^^Z 



0.43+iln(^^/Te) 



-F0.469(£;^/Te)-'/^ 



cm / sec. 
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For 1 eV < Te/Z^ < 15 eV, this becomes approximately"^^ 

(14) ariZ) = 2.7 X 10"^^Z^Te"^/^ cm^/sec. 

Collisional (three-body) recombination rate for singly ionized plasma 

(15) as = 8.75 x 10"^"^ Tg"^'® cm® /sec. 

Photoionization cross section for ions in level n, I (short- wavelength limit): 

(16) cTph(n, = 1.64 X 10"^®Z^/n^K^+^^ cm^, 

where K is the wavenumber in Rydbergs (1 Rydberg = 1.0974 x 10^ cm~^). 



55 



Ionization Equilibrium Models 

Saha equilibrium:^^ 



(17) — \/ ^ \ = 6.0 X 10^^ ^ , exp ' cm ^ 

iV,*(Z-l) ^f-i ' 



where g'^ is the statistical weight for level n of charge state Z and (n 
is the ionization energy of the neutral atom initially in level {n,l), given 
Eq. (2). 

In a steady state at high electron density, 

NeN*{Z) S{Z - 1) 



(18) 



N*(Z - 1) as 



a function only of T. 

Conditions for LTE:^® 

(a) Collisional and radiative excitation rates for a level n must satisfy 

(19) Yr,m Z 10Ar^m. 

(b) Electron density must satisfy 

(20) Ne>7x 10^^Z^n-^^/^(T/£;2>)^^^cm-^ 
Steady state condition in corona model: 

N(Z-l) Oir 



(21) 



N{Z) S{Z-1) 
Corona model is applicable if*° 

(22) lO^^tj"^ < Ne < 10^®Te^/^ cm"^ 

where tj is the ionization time. 
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Radiation 

N. B. Energies and temperatures are in eV; all other quantities are in 
cgs units except where noted. Z is the charge state (Z = refers to a neutral 
atom); the subscript e labels electrons. N is number density. 

Average radiative decay rate of a state with principal quantum number n is 



(23) A„ = ^ 



Anm = 1-6 X 10^°Z^n"^/^sec. 

m < n 



Natural linewidth {AE in eV): 

(24) A£; At = h = 4.14 X 10"^^ eVsec, 

where At is the lifetime of the line. 
Doppler width: 

(25) AA/A = 7.7 X 10"^(T///)^/^, 

where /j, is the mass of the emitting atom or ion scaled by the proton mass. 
Optical depth for a Doppler-broadened line:^^ 

(26) T = 3.52 X 10~^^ Um MM /kT)^^'^ NL = 5.4 x 10~^ fmnHl^/T)^^^ NL, 

where fnm is the absorption oscillator strength, A is the wavelength, and L is 
the physical depth of the plasma; M , N , and T are the mass, number density, 
and temperature of the absorber; ji is M divided by the proton mass. Optically 
thin means r < 1. 

Resonance absorption cross section at center of line: 

(27) (^\=\c = 5-6 X 10"^^A^/AAcm^. 

Wien displacement law (wavelength of maximum black-body emission): 

(28) Amax = 2.50 X 10"^T"^ cm. 

Radiation from the surface of a black body at temperature T: 

(29) W ^ 1.03 X lO'^T^ watt/cm^. 
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Bremsstrahlung from hydrogen- like plasma: 

(30) Psr = 1.69 X 10"^^A/"eTe^/^ ^ [z^Ar(Z)] watt/cm^, 

where the sum is over all ionization states Z. 
Bremsstrahlung optical depth:'*^ 

(31) r = 5.0 X 10~^^ NeNiZ^gLT'^^"^ , 

where 'g ^ 1.2 is an average Gaunt factor and L is the physical path length. 

Inverse bremsstrahlung absorption coefficient'*^ for radiation of angular fre- 
quency uj: 

(32) K = 3.1 X 10"^Zne^ In A T"^/^cj"^ (1 - cj^/cj^)"^/^ cm"^ 

here A is the electron thermal velocity divided by V, where V is the larger of 
u! and ujp multiplied by the larger of Ze^/kT and h/^rnkT)^^"^. 

Recombination (free-bound) radiation: 



(33) Pr = 1.69 X lO'^^ATeTe^/^ ^ Z^Ar(Z) 



oo 

Te J 



watt / cm" 
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.41 



Cyclotron radiation^® in magnetic field B: 

(34) Pc = 6.21 X lO'^^B^ATeTe watt/cm^. 

For ATefcTe = NikTi = S^/167r (/3 = 1, isothermal plasma), 

(35) Pc = 5.00 X 1Q~^^nItI watt/cm^. 

Cyclotron radiation energy loss e-folding time for a single electron 

, , 9.0 X lO^S-2 

(36) tc ^ -— sec, 

2.5 + 7 

where 7 is the kinetic plus rest energy divided by the rest energy mc . 
Number of cyclotron harmonics^^ trapped in a medium of finite depth L: 

(37) mtr = (57/3BL)'/^ 
where P = SnNkT/B"^. 

Line radiation is given by summing Eq. (9) over all species in the plasma. 
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ATOMIC SPECTROSCOPY 



Spectroscopic notation combines observational and theoretical elements. 
Observationally, spectral lines are grouped in series with line spacings which 
decrease toward the series limit. Every line can be related theoretically to a 
transition between two atomic states, each identified by its quantum numbers. 

Ionization levels are indicated by roman numerals. Thus C I is unionized 
carbon, C II is singly ionized, etc. The state of a one-electron atom (hydrogen) 
or ion (Hell, Li III, etc.) is specified by identifying the principal quantum 
number n — 1, 2, ... , the orbital angular momentum I — 0, 1, . . . , n — 1, and 
the spin angular momentum s = ±i. The total angular momentum j is the 

magnitude of the vector sum of 1 and s, j' = i ± ^ (j > ^). The letters s, 
p, d, f, g, h, i, k, 1, . . . , respectively, are associated with angular momenta 
I = 0, 1, 2, 3, 4, 5, 6, 7, 8, .... The atomic states of hydrogen and hydrogenic 
ions are degenerate: neglecting fine structure, their energies depend only on n 
according to 

_ RoohcZ^n~^ _ RyZ^ 
~ 1 + m/M " n2 ' 

where h is Planck's constant, c is the velocity of light, m is the electron mass, 
M and Z are the mass and charge state of the nucleus, and 

Roo = 109, 737cm"^ 

is the Rydberg constant. If En is divided by he, the result is in wavenumber 
units. The energy associated with a transition m — >■ n is given by 

AEmn = Ry(l/m^ - 1/n^), 

with m < n (m > n) for absorption (emission) lines. 

For hydrogen and hydrogenic ions the series of lines belonging to the 
transitions m — >■ n have conventional names: 



Transition 


1— )-n 2— )-n 3— )-n 4— )-n 5— )-n 6— )-n 


Name 


Lyman Balmer Paschen Brackett Pfund Humphreys 



Successive lines in any series are denoted a, /3, 7, etc. Thus the transition 1 — >■ 
3 gives rise to the Lyman-/3 line. Relativistic effects, quantum electrodynamic 
effects (e.g., the Lamb shift), and interactions between the nuclear magnetic 
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moment and the magnetic field due to the electron produce small shifts and 
splittings, ^ 10~^cm~^; these last are called "hyperfine structure." 

In many-electron atoms the electrons are grouped in closed and open 
shells, with spectroscopic properties determined mainly by the outer shell. 
Shell energies depend primarily on n; the shells corresponding to n = 1, 2, 
3, . . . are called K , L, M , etc. A shell is made up of subshells of different 
angular momenta, each labeled according to the values of n, I, and the number 
of electrons it contains out of the maximum possible number, 2(21 + 1). For 
example, 2p^ indicates that there are 5 electrons in the subshell corresponding 
to Z = 1 (denoted by p) and n — 2. 

In the lighter elements the electrons fill up subshells within each shell 
in the order s, p, d, etc., and no shell acquires electrons until the lower shells 
are full. In the heavier elements this rule does not always hold. But if a 
particular subshell is filled in a noble gas, then the same subshell is filled in 
the atoms of all elements that come later in the periodic table. The ground 
state configurations of the noble gases are as follows: 

He Is^ 

Ne ls^2s^2p^ 

Ar ls^2s^2p^3s^3p^ 

Kr Is22s22p63s23p^3d^°4s24p6 

Xe Is22s22p63s23p63d^04s24p64d^05s25p^ 

Rn Is22s22p63s23p63di04s24p64di04fi45s25p65di°6s26p6 

Alkali metals (Li, Na, K, etc.) resemble hydrogen; their transitions are de- 
scribed by giving n and I in the initial and final states for the single outer 
(valence) electron. 

For general transitions in most atoms the atomic states are specified in 
terms of the parity ( — 1)^'^ and the magnitudes of the orbital angular momen- 
tum L = El^, the spin S = Ss^, and the total angular momentum J = L + S, 
where all sums are carried out over the unfilled subshells (the filled ones sum 
to zero). If a magnetic field is present the projections Ml, Ms, and M of 
L, S, and J along the field are also needed. The quantum numbers satisfy 
\Ml\ < L < vl, \Ms\ < S < v/2, and |M| < J < L + 5, where v is the 
number of electrons in the unfilled subshell. Upper-case letters S, P, D, etc., 
stand for L = 0, 1, 2, etc., in analogy with the notation for a single electron. 
For example, the ground state of CI is described by 3p^ ^^3/2" The first part 
indicates that there are 5 electrons in the subshell corresponding to n = 3 and 
I = 1. (The closed inner subshells ls^2s^2p^3s^, identical with the configura- 
tion of Mg, are usually omitted.) The symbol 'P' indicates that the angular 
momenta of the outer electrons combine to give L = 1. The prefix '2' repre- 
sents the value of the multiplicity 2S -\- 1 (the number of states with nearly the 
same energy), which is equivalent to specifying S = ^. The subscript 3/2 is 
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the value of J. The superscript 'o' indicates that the state has odd parity; it 
would be omitted if the state were even. 

The notation for excited states is similar. For example, helium has a state 
ls2s ^Si which lies 19.72eV (159, 856 cm~^) above the ground state Is^ ^Sq. 

But the two "terms" do not "combine" (transitions between them do not occur) 
because this would violate, e.g., the quantum-mechanical selection rule that 
the parity must change from odd to even or from even to odd. For electric 
dipole transitions (the only ones possible in the long- wavelength limit), other 
selection rules are that the value of I of only one electron can change, and only 
by AZ = ±1; AS' = 0; AL = ±1 or 0; and AJ = ±1 or (but L = does not 
combine with L = and J = does not combine with J = 0). Transitions 
are possible between the helium ground state (which has = 0, L = 0, J = 0, 
and even parity) and, e.g., the state ls2p ""^P^ (with 5 = 0, L = l, J = l, 
odd parity, excitation energy 21.22 eV). These rules hold accurately only for 
light atoms in the absence of strong electric or magnetic fields. Transitions 
that obey the selection rules are called "allowed" ; those that do not are called 
"forbidden." 

The amount of information needed to adequately characterize a state in- 

A O 

creases with the number of electrons; this is reflected in the notation. Thus 
O II has an allowed transition between the states 2p 3p 
^■^7/2 2p^("'^D)3d'' ^F7/2 (and between the states obtained by changing 

J from 7/2 to 5/2 in either or both terms). Here both states have two elec- 
trons with n = 2 and / = 1; the closed subshells ls^2s^ are not shown. The 
outer (n = 3) electron has Z = 1 in the flrst state and Z = 2 in the second. 
The prime indicates that if the outermost electron were removed by ionization, 
the resulting ion would not be in its lowest energy state. The expression (""^D) 
give the multiplicity and total angular momentum of the "parent" term, i.e., 
the subshell immediately below the valence subshell; this is understood to be 
the same in both states. (Grandparents, etc., sometimes have to be specified 
in heavier atoms and ions.) Another example"^^ is the allowed transition from 
2p2(3p)3p (or ^Pg/s) to 2p2(^D)3d' ^Si/2, in which there is a "spin 

flip" (from antiparallel to parallel) in the n = 2, Z = 1 subshell, as well as 
changes from one state to the other in the value of I for the valence electron 
and in L. 

The description of flne structure. Stark and Zeeman effects, spectra of 
highly ionized or heavy atoms, etc., is more complicated. The most important 
difference between optical and X-ray spectra is that the latter involve energy 
changes of the inner electrons rather than the outer ones; often several electrons 
participate. 
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COMPLEX (DUSTY) PLASMAS 



Complex (dusty) plasmas (CDPs) may be regarded as a new and unusual 
state of matter. CDPs contain charged microparticles (dust grains) in addition 
to electrons, ions, and neutral gas. Electrostatic coupling between the grains 
can vary over a wide range, so that the states of CDPs can change from weakly 
coupled (gaseous) to crystalline. CDPs can be investigated at the kinetic level 
(individual particles are easily visualized and relevant time scales are accessi- 
ble). CDPs are of interest as a non-Hamiltonian system of interacting particles 
and as a means to study generic fundamental physics of self-organization, pat- 
tern formation, phase transitions, and scaling. Their discovery has therefore 
opened new ways of precision investigations in many-particle physics. 



Typical experimental dust properties 

grain size (radius) a c^i 0.3 — 30 ^tm, mass md 3 X 10~^ — 3 X lO""*^^ g, number 
density (in terms of the interparticle distance) ~ 
temperature ~ 3 X 10~^ - 10^ eV. 

Typical discharge (bulk) plasmas 



10^ - 10^ cm"^. 



gas pressure p ~ 10 — 1 Torr, Tj ~ T„ ~ 



3 X 10~^ eV, ~ 7 X 10* cm/s 
(Ar), Te ~ 0.3 - 3 eV, ~ rie ~ 10^ - 10^° cm~^, screening length Ad ~ 
Az5i ~ 20 - 200 /xm, u;^^ ~ 2 x 10^ - 2 x 10^ s~^ (Ar). B fields up to S ~ 3 T. 



Dimensiotiless 

Havnes parameter 

normalized charge 

dust-dust scattering parameter 

dust-plasma scattering parameter 

coupling parameter 

lattice parameter 

particle parameter 

lattice magnetization parameter 

Typical experimental values: P lO"'*- 
charges), P < 10^, k ~ 0.3 - 10, a ~ 10" 

Frequencies 

dust plasma frequency 

charge fluctuation frequency 



P = \Z\nd/ne 
z = |Z|e^/fcTea 
= Z^e'^/kTdXD 

Pe,i = \Z\e'^ /kTe,i\D 

r = (Z2eVfcTdA)exp(-A/Ai3) 
K = A/Ad 
a = a/A 
M = A/rd 



-10^,z ~ 

-4 



2-4 (Z ^ 10^-10^ electron 



3 X 10 



oJpd = (47rZ^e^nd/md)^/^ 
- (\Z\^mi/md)^^^uj 



pi 



1 + z 
V2^ 



{a/XD)oj 



pi 
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dust-gas friction rate 
dust gyrofrequency 

Velocities 

dust thermal velocity 

dust acoustic wave velocity 

dust Alfven wave velocity 
dust-acoustic Mach number 
dust magnetic Mach number 

dust lattice (acoustic) wave velocity 



h>nd ~ lOa^p/mdVTr 
LOcd = ZeB/rudC 



CdA = (^pd^D 

- (\Z\j^mi/mdy/'^VTi 
VAd = -B/(47rndmd)^/^ 
V/Cda 

V/vAd 



The range of the dust-lattice wavenumbers is KA < tt The functions Fi^t^K.) 
for longitudinal and transverse waves can be approximated^'*''*^ with accuracy 
< 1% in the range k < 5: 

Fi ~ 2.70k^/^(1 - 0.096k - 0.004k^), Ft 0.51k(1 - 0.039k^), 



Lengths 

frictional dissipation length 

dust Coulomb radius 
dust gyroradius 

Grain Charging 



Lu = VTd/^nd 

Rce,i = \Z\e'^/kTe 

Vd = VT^/(^cd 



The charge evolution equation is d\Z\/dt = li — le- From orbital motion 
limited (OML) theory^® in the collisionless limit len(in) ^ ^ o,: 



/g = v87ra UeVTe exp( — z), li = v87ra riiVT- I 1 -\ z 



Grains are charged negatively. The grain charge can vary in response to spatial 
and temporal variations of the plasma. Charge fluctuations are always present, 
with frequency cuch- Other charging mechanisms are photoemission, secondary 
emission, thermionic emission, field emission, etc. Charged dust grains change 
the plasma composition, keeping quasineutrality. A measure of this is the 
Havnes parameter P — \Z\nd/ne- The balance of 1^ and li yields 



exp( — 2;) = 



rrii T, 



1/2 



T, 



T, 
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When the relative charge density of dust is large, P ^ 1, the grain charge Z 
monotonically decreases. 

Forces and momentum transfer 

In addition to the usual electromagnetic forces, grains in complex plasmas are 
also subject to: gravity force Fg = rridS', thermophoretic force 

4^277 , 2 / 

Fth 7^{(^ /vTn)l^nVTn 

15 

(where Kn is the coefficient of gas thermal conductivity); forces associated 
with the momentum transfer from other species, Fq, = —mdi^adL^ ad-, i-e., 
neutral, ion, and electron drag. For collisions between charged particles, two 
limiting cases are distinguished by the magnitude of the scattering parameter 
I^Q,. When /3q, <C 1 the result is independent of the sign of the potential. When 
/3q, ^ 1, the results for repulsive and attractive interaction potentials are 
different. For typical complex plasmas the hierarchy of scattering parameters 
is /3e(-- 0.01 - 0.3) < 1-30) <C I3d{^ 10^ - 3 x 10^). The generic 

expressions for different types of collisions are^^ 

i^ad = (4v27r/3)(mc/md)a nocVToc^ad 



Electron-dust collisions 



1 2 

<l>ed ^ -Z Aed /3e < 1 



Ion-dust collisions 



^ ^ f \z'^{T^lTifA,d /3i < 5 

\2(AD/a)='(ln2/3i + 21n/3i+2), /3i > 13 

Dust-dust collisons 

= [ ^d^'^d fid « 1 

\(Az5/a)2[ln4/3d-lnln4/3d], /3d > 1 
where Zd = Z e jakTd- 

For Vdd ^nd the complex plasma is in a two-phase state, and for u^d ^ ^dd 
we have merely tracer particles (dust-neutral gas interaction dominates). The 
momentum transfer cross section is proportional to the Coulomb logarithm 
Actd when the Coulomb scattering theory is applicable. It is determined by 
integration over the impact parameters, from pmin to pmax- Pmin is due to finite 
grain size and is given by OML theory, pmax = for repulsive interaction 

(applicable for /3a <C 1), and pmax = Ai3(l + 2;0q:)^^^ for attractive interaction 
(applicable up to /3q, < 5). 
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For repulsive interaction (electron- dust and dust-dust) 




where z^. = z, = a, and ad = 2a. 
For ion-dust (attraction) 




l + 2(T,/Te)(Azj/a)x 
l + 2{Ti/Te)x 



ln(2x — l)dx, 



dx. 



For Vdd ^ i^nd the complex plasma behaves like a one phase system (dust-dust 
interaction dominates). 

Phase Diagram of Complex Plasmas 

The figure below represents different "phase states" of CDPs as functions of 
the electrostatic coupling parameter P and k or a, respectively. The verti- 
cal dashed line at k = 1 conditionally divides the system into Coulomb and 
Yukawa parts. With respect to the usual plasma phase, in the diagram be- 
low the complex plasmas are "located" mostly in the strong coupling regime 
(equivalent to the top left corner). 

Regions I (V) represent Coulomb (Yukawa) crystals, the crystallization condi- 
tion is^^ F > 106(1 -I- K -I- fv^/2)~^. Regions II (VI) are for Coulomb (Yukawa) 
non-ideal plasmas — the characteristic range of dust-dust interaction (in terms 
of the momentum transfer) is larger than the intergrain distance (in terms of 

the Wigner-Seitz radius), (ct/tt)-'-/^ > (47r/3)~-'-/^A, which implies that the 
interaction is essentially multiparticle. 



a 



=A/a 



10' 



10^ 



10 



10' 



10^ 



U 10": 
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Regions III (VII and VIII) correspond to 
Coulomb (Yukawa) ideal gases. The range 
of dust-dust interaction is smaller than the 
intergrain distance and only pair collisions 
are important. In addition, in the region 
VIII the pair Yukawa interaction asymp- 
totically reduces to the hard sphere limit, 
forming a "Yukawa granular medium". In 
region IV the electrostatic interaction is 
unimportant and the system is like a uaual 
granular medium. 
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